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ABSTRACT 


Solid/liquid distribution coefficients (weight basis) have been experimentally 

determined for a number of trace elements for olivine, orthopyroxene, plagioclase 

and ilmenite. Values of di^fibution coefficients measured at 1200°C and f_ of 

O2 

-13.0 

10 for liquids similar in composition to the olivine-opx-plagioclase peritectic 

in the pseudoternary system (Fe,Mg) 2 Si 0 ^-CaAl 2 Si 20 g-Si 02 are as follows: 

olivine; Ce = 0.010 ±.007, Sm = 0.015 ±.005, Eu = 0.015 ±.005, Yb = 0.033 ±.002, 

Cr = 1.2 ±.l; orthopyroxene; Ba = 0,011 ±.005, Cr = 5.2 ±.2: anorthite; Rb = 0.017 ±.008, 

Ba = 0.15 ±.03. Values measured at 1140®C and f„ of 10 for liquids similar 

O2 

in composition to high-Ti mare basalts are as follows: Ilmenite; Ce = 0.006 ±.003, 

Sm = 0.010 ±.002, Eu = 0.007 ±.003, Yb = 0.075 ±.005. Cr distribution coefficients 
decrease with decreasing oxygen fugacity for both pyroxene and olivine, but are 
larger for pyroxene than for olivine at all fugacities investigated. The variation 
of D with oxygen fugacity indicates that a substantial fraction of the Cr 
is divalent at lunar oxygen fugacities. 

Major and trace element partitioning and relevant phase equilibria are used 
to investigate possible parent-daughter relationships between a number of highland 
samples and highly evolved KREEP-rich materials. Out of about 80 highland samples 
tested, 33 are found to be possible parents to the KREEP-rich materials. The 
average composition of these 33 samples is very similar to that of the Low-K Era 
Mauro basalt (LKFM) . A model is proposed in which LKEM-type material was produced 
by fractionation of large amounts of olivine and lesser amounts of plagioclase 
(and possibly other minor phases) from undifferentiated lunar material at 4.4 M. 

During the same episode of igneous activity, some LKEM material differentiated to 
form a series of more evolved members of the KREEP suite. Including material of 
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peritectic bulk composition in the (Mg,Fe) 2 Si 0 ^-CaAl 2 Si 20 g-Si 02 pseudoternary 
system. In a later igneous event at 3.9 jSE, some of this peritectic material 
further differentiated, producing more highly evolved KREEP such as sample 15386. 

The model is consistent with phase equilibria, major and trace element partitioning, 
and Rb-Sr isotopic data. 
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INTRODUCTION 


There is a large family of lunar highland rocks with marked enrichment in a 
number of minor and trace elements (K, P and the rare earth elements, among others) 
which are not readily incorporated into the crystal lattices of the major 
rock-forming minerals. These KREEP-rich highland rocks are also characterized by 
bulk compositions which, when plotted in the (Mg,Fe) 2 Si 0 ^-Si 02 -CaAl 2 Si 20 g 
pseudoternary liquidus diagram, form a continuum extending along the olivine- 
plagioclase cotectic line through the olivine-plagioclase-pyroxene peritectic point 
and along the plagioclase-pyroxene cotectic line (Fig. 1 and Table 1). A few of 
these rocks, such as 15386, have igneous textures. Most of them are breccias, and 
many of these are polymict. Most KREEP-rich samples are therefore mechanical mixtures 
of fragments of rocks or minerals. As such, they could be considered mixtures of 
the “essence of KREEP or the "KREEP component" with other lunar rock types such 
that the high overall concentration and the relative abundances of K, REE, P 
and associated "incompatible" elements are dominated by the KREEP ■ component, while 
the bulk composition is determined by the relative proportions of all the lunar 
rock types in the mixture. According to this point of view differences in the major 
element composition of KREEP-rich samples are primarily determined by a variable 
history of mechanical mixing. 

Other interpretations are possible. Walker ^ (1972) pointed out that the 

major element compositions of the KREEP-rich rocks cluster along solid-liquid 
saturation curves of the pseudoternary phase diagram, suggesting that they had 
evolved as a result of liquid separations from equilibrated partially molten phase 
assemblages within the pseudoternary system. In approaching the problem of KREEP 
we assume that the association of their bulk chemistry with liquid saturation 
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boundaries is not coincidental. Therefore, even though these rocks have been 
subjected to impact fragmentation, recombination and thermal metamorphic effects 
on the lunar surface we take the view that their bulk compositions have not been 
greatly altered and that the general characteristics of their chemical compositions 
can be interpreted mainly in terms of equilibrium crystallization and melting. 

It is consistent with the major element compositions of all these rocks that 
they are related by the solid-liquid phase equilibria depicted in Fig. 1. Therefore, 
it is logical to ask whether such relationships are also compatible with the trace 
element chemical data. Using the phase equilibria in Fig. 1 and the solid/liquid 
distribution coefficients for a nxomber of key elements it is possible to relate the 
proportions of solid and liquid during differentiation to the resulting trace 
element fractionation trends. Conversely, given the trace element concentrations 
in an assumed parent assemblage and its derivative, we can calculate the proportions 
of the major phases involved in the fractionation process and check for compatibility • 
with the phase diagram. Using the above approach Weill and McKay (1975) demonstrated 
that -it is improbable that the KEEEP-rich rocks were derived, by igneous fractionation 
from what is now estimated to be .the average highlands composition (Taylor and Jakes, 
1974). Weill and McKay (1975) also investigated the possibility of deriving KREEP 

r 

liquids by one-step equilibrium partial melting or crystallization of compositions 
similar to those proposed for the whole moon (Ganapathy and Anders, 1974; Taylor 
and Jakes, 1974) . It was concluded that such an origin is consistent with the 
known constraints of phase equilibria and trace element distribution behavior if 
a few percent or less of KREEP liquid^equilibrated with a mineral assemblage of 
olivine, Ca-poor pyroxene and plagioclase. 

We believe these conclusions are still valid, but the resulting model for the 
origin of KSEEP-rich rocks is not unique or complete. Our previous calculations 
addressed themselves to the derivation of KREEP-rich compositions corresponding 
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closely to the pseudoternary peritectic and did not attempt to investigate the 
possible relationships among the total range of KREEP-rich rocks. Perhaps a 
multiple-step origin is equally plausible (e.g. , recent Sm/Nd isotopic data 
presented by Lugmair et al. (1976) indicate a minimum of two steps for the origin 
of KEIEEP sample 15382). Hence, it is appropriate in this paper to consider the 
possibility of multiple-step derivation of KREEP and the relationships between 
the KREEP-rich compositions generated at successive stages of evolution. In this 
context we also discuss the effect of the igneous processes involved in the 
evolution of KREEP on Rb-Sr isotopes. Finally, we also present some recent 

experimental determinations of solid/liquid distribution coefficients for the REE 
in olivine and ilmenite and Cr, Ba and Rb in plagioclase and pyroxene. 


EXPERIMEHTAL METHODS 

A set of bulk compositions were prepared so that when they were held at 1200°C 
and at oxygen fugacities low enough to maintain the Fe predominantly in the 
divalent state they produced several percent of either olivine, orthopyroxene or 
plagioclase in equilibrium with liquids ofu compositions close to the peritectic 
point in Fig. 1, These bulk compositions were then enriched with an additional 
1-3 wt % of one of the large-ion-lithophile elements (or about 0.5 wt Z in the case 
of Cr) and held at temperature and oxygen fugacity (H 2 /CO 2 mixing) in SiC resistance 
furnaces for 50—150 hours before quenching. Ilmenite/liquid distribution coefficients 
were determined in bulk systems corresponding to high-Ti mare basalts with 1-3 wt % 
LIL element added. In order to minimize Fe losses during runs the charges were held 
in PtRh wire loops which had previously been held in contact with similar 
compositions at 

crystals and glass were analyzed by electron probe. The distribution coefficients 
and the details of phase compositions for these experiments are summarized in 


T and f_ similar to run conditions. After quenching, 
O2 


Table 2 
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EXPERIMENTAL RESULTS 


Values o£ the solid/liquid distribution coefficients determined in this study- 

are given in Table 2 , along with some previously determined coefficients from Weill 

and McKay (1975) , These distribution coefficients have been determined specifically 

for the range of composition and temperature appropriate to equilibria along the 

plagioclas e-olivine cotectic and near the plagioclase-olivine-pyroxene peritectic in 

the system (Mg,Fe) SiO -CaAl Si 0 -SiO • Attempts to apply the coefficients to systems 
2. H A A Q A 


with substantially different phase chemistry involves an uncertain extrapolation 

which should be approached with caution. Only the plagioclase/liquid distribution 
' ■ orthopyroxene 

coefficients for Eu and the /liquid and olivine /liquid coefficients for Cr show 

Cr 

significant dependence on oxygen fugacity. In Eig. 2 we have plotted 

D^l/liq as a function of oxygen fugacity. As in the case of Eu in plagioclase and 

Ca-pyroxene (Weill and Drake, 1973; Weill et al. 1974), these data can be most 


easily interpreted in terms of variations of the fractions of the cation, in this 
case Cr, in the divalent and trivalent states. These variations may be expressed as 


+2 +3 

an overall Cr distribution coefficient, D = D 2 (Cr ) + ^^C^r ). According to this 

interpretation the data plotted in Fig. 2 indicate that > D^ both for olivine and 

orthopyroxene. The difference is more pronounced in orthopyroxene. Schreiber and 

Haskin (1976) have published some experimental distribution coefficients for Cr in 

enstatite?. dlopside and forsterite. Their determinations were carried out using 

iron-free synthetic compositions including one (their composition FAS) lying within 

the ternary system Mg SiO.-CaAl Si„0_-Si0 . Composition FAS has some obvious 

Ah a a a 


similarities to the compositions used in this study (near the peritectic point 

in the pseudoternary system (Mg,Fe) „SiO .-CaAl„Si_On-SiO_) although compositions 

2 4 2 2 o 2 


in the iron-free. system have higher liquidus temperatures than the corresponding Fe- 
bearing liquids lying at the same point in the pseudoternary diagram. Because of the 
differences between the compositions investigated by Schreiber and Haskin and 
those which we investigated, and because the lowest temperature for which they 
determined Cr distribution coefficients (1300 “C) is considerably higher than the 
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temperature at which we determined coefficients (1200'’C) , their results cannot he 

directly compared with ours. Nevertheless, the results of the present study have 

several Important general features in common with the results of Schreiher and 

Oir 

Haskin for their composition FAS: both studies show D decreasing with decreasing 

f^^ < 10 in both orthopyroxene and olivine, with the magnitude 

of the decrease being much larger for orthopyroxene; we infer from the variation 

of the D's with f^ that a substantial fraction of the Cr in our experimental liquids 

^ -13 

is divalent at 1200“C and f of 10 (we are unable to quantitatively determine 

2 

the value of this fraction) , while Schreiher and Haskin find that more than 75% 


of the Cr in their liquids (as measured by titration and EPR) is divalent at 1300°C 
-10 

and f^^ of 10 ; we find that for orthopyroxene > 9 at 1200®C, while Schreiher 

Ftt r*T" 

and Haskin find that D- > 7 at 1300“C. Values of D , , D , „Cr 

3 pl/liq’ px/lxq and 


appropriate to the lunar range of T-f^ are given in Table 2. 

°2 

DISCUSSION 

A highly evolved KREEP sample - 25386 

In terms of major element composition, sample 15386 is the most evolved member 

of the KREEP suite shown in Fig. 1. Its bulk composition places it along the 

plagioclase-pyroxene saturation curve, down-temperature from the olivine-pyroxene- 

plagioclase peritectic point. It is one of the relatively few samples of KREEP 

with unmistakably igneous texture, and its low Ni content (Nyquist et al, 1975), 

Implying a low amount of meteoritic contamination, further enhances the likelihood 

that its immediate origins were igneous. According to the phase equilibrium diagram 

of Fig. 1, a liquid of composition 15386 is an intermediate product of differentiation 

(via either partial melting or crystallization) of a parental composition near the 

peritectic point. The pseudoternary phase diagram and the distribution coefficient 
Fe/Mg 

%x/liq combined with the known Mg and Fe contents of 15386 to determine the 

major element concentrations in a model parent P (Table 1 and Fig. 5) and the relative 
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praportions of coexisting plagioclase, pyroxene and liquid 15386. Fortunately 
15386 is a very well documented Apollo 15 crystalline KREEP sample. In addition to 
its major element content, there are reliable trace element and whole rock Rb/Sr 
data, as well as an internal Rb/Sr isochron age. The known trace element concentrations 
in 15386 allow us to calculate the trace element concentrations in parental composition 
P according to the mass balance equation for each trace element: 

D - Z D.W. = (C - C^F)/(1 - F)a [1] 

.-XI O L Ij 

X=1 


where D is the composite solid/liquid distribution coefficient, is the 

distribution coefficient and W, is the weight fraction of mineral i 
n ^ 

= 1) , while F is the weight fraction of the equilibrated liquid relative to 
the total weight of the system. C is the concentration of the trace element in the 
liquid (i.e. 15386) and is the concentration of the trace element in the total 


system (i.e,, parental composition P) . In Table 1 we have listed the major and trace 


element data for 15386 and its model parent composition P. Plagioclase is involved 
in the differentiation step P -s- 15386, and because very nearly one 

hundred times greater than » this step also involves an increase of Rb/Sr 

from 0.061 -in^composition P .to 0,099 in 15386. We will return to this point in 
discussing the Rb/Sr isotopes. 


PossihZe Tetati-onship hetween pev-ttectic KREEP and othev highland rocks 

The peritectic point of the pseudoternary phase diagram is obviously an 
important composition. Many of the KREEP rocks plot close to this point, and we 
have also just discussed the possibility that such liquids are parental to some 
of the more evolved KREEP samples such as 15386. Mineral assemblages of olivine + 
plagioclase + pyroxene produce initial melts of the peritectic composition. 

Conversely, during crystallization in the "highlands" systems residual liquid 
compositions approach the peritectic. It is desirable therefore to consider the question 
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of the possible origin of KB.EEP peritectic liquids and examine whether there are 
any rock compositions in the highlands which could have given rise to KREEP-rich 
peritectic liquids via either equilibrium crystallization or partial melting. In 
pursuing this question we will use two typical examples of KREEP ’’peritectic” 
compositions; the possible parent of 15386 which we have just calculated (model 
composition P) and the soil 14163 which is representative of Apollo 14 ICREEP 
breccias and plots very close to the peritectic in Fig. 1. Both of these com- 
positions are given in Table 1. In keeping with our previous discussion we assume 
that the evolution to liquids P and 14163 was governed by liquid-solid differentiation 
and' that the trace element partitioning is described by the mass balance of equation 
[1] . The refractory solid phase assemblages involved are combinations of plagioclase, 
pyroxene and olivine. The peritectic liquid is assumed to be essentially equilibrated 
with these mineral assemblages either after partial melting or crystallization of the 
parent material. As candidates for parental materials to peritectic liquids we have 
considered highland rocks (approximately 80) for which we could find published major 
and trace element analyses of the same subsample. 

There are many sets and subsets of abundance criteria which could be used as 
filters to select the potential- parents of KREEP-rich 

peritectic liquids. Our distribution coefficient data (Table 2) Include 
typical "excluded" (i.e., trivalent) REE, an atypical REE (Eu) and the important 
major elements Mg and Fe. The approach we have employed is first to determine, 
for each of the parental candidates, the relative -proportions of each phase which 
must be present in the equilibrated assemblage (i.e., peritectic liquid plus solids) 
in order that the liquid have the same content of the excluded trace element 8m 
as that observed in 14163 or composition P. We then check whether these relative 
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phase proportions would also yield a liquid with the abundances of Eu and 
Mg/ (Mg + Fe) observed in 14163 or composition P. The calculations proceed as follows: 
(1) We first assign a value of D = 0 for Sm. This allows an approximate calculation 
of F according to equation 1. At this point we arbitrarily reject candidates for 
which F > 0.8 since such materials are already so enriched in KREEP components that 
they are of little interest in considering the parentage of KREEP. For the 
remaining candidates, we use the approximate value of F to calculate the approximate 

bulk composition (and hence the mode) of the residual solid assemblage according to 
the simple mass balance relationship, = FG^^ + (1 - F)Cg , where 

is the observed concentration of a major element in the potential parent, C, 

Li 

is the concentration in 14163 or P, and C_ is the concentration in the 

D 

residual solid assemblage. The approximate mode in turn allows a refined value of 

Dg^ to be calculated (equation 1) using the individual mineral distribution 

coefficients for Sm from Table 2. The refined is in turn used to calculate a 

Sm 

refined value for F, etc., etc. The iterations are terminated when successive solutions 
of F differ by less than 10% relative. (2) The final value for _F_and the corresponding 
mode of the refractory solid assemblage are then combined with the distribution 
coefficients for Eu, Mg and Fe (Table 2) in equation [1] to calculate the Eu and 
Mg/ (Mg + Fe) abundances in the liquid. 

The calculated abundances of Eu and Mg/ (Mg + Fe) for each candidate parent 

are plotted in Fig. 3 where they are compared to 14163. The rectangle around 14163 

0 

represents estimates of the combined uncertainties in the distriibution coefficients 
of Sm, Eu, Fe and Mg. Highland rocks which are potential parental compositions 
for 14163 must plot within the rectangle of Fig. 3. Note that the average highland 
composition of Taylor and Jakes (1974) produces a liquid with the appropriate 
Sm and Eu content but too little Mg/ (Mg + Fe) as pointed out by Weill and McKay (1975) . 
We next consider whether there is anything distinctive about the bulk compositions 
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of samples which qualify as potential parents to 14163. In Fig. 4a we have plotted 
all candidates tested. The compositions span nearly the entire range of highland rock 
types. In Fig 4h we have plotted all candidates that meet the criteria for parental 
material. The large majority of them are close to the olivine-plagioclase cotectic. 

A similar set of calculations was repeated using composition P as the daughter liquid 
instead of 14163, and the results were very similar. Since the criteria were com- 
pletely independent of the location of the phase saturation boundaries, the results 
may be Interpreted as independently reinforcing the notion that many of the KREEP suite 
of rocks are related by the phase equilibria depicted in Fig. 1. 

In all, 33 highland rock compositions were found to be potential parents 

for peritectic liquids P or 14163. The average composition of all these potential < 
parents is given as composition C in Table 1. This average is very similar to the 
Apollo 15 Low-K Fra Mauro (LKFM) basalt type (also shown in Table 1) proposed 
by Reid (1972) . Taylor and Jakes (1974) concluded on the basis of orbiter 

data that the LKFM composition is an abundant constituent of the lunar highlands. 

Thus, the constraints of low pressure phase equilibria and of major and trace element 
partitioning are satisfied by models which relate KREEP peritectic liquids 
(and their possible derivatives such as 15386) to an abundant crustal 
rock type via equilibrium crystallization and partial melting. 

Fossibte ovigin of KREEP 

We have shown that it is possible to consider the KREEP-rich rocks as an igneous 
suite evolved by solid-liquid differentiation along the crystal saturation boundaries 
of the highlands pseudoternary phase equilibrium diagram. The least evolved 
members of this suite plot along the olivine-plagioclase cotectic with an average 
composition corresponding to point C in Fig. 5. The ultimate origin of such, 
liquids is a question which may never be answered and certainly will be extremely 
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difficult to answer in any detail. Nevertheless it is tempting to consider that 
a liquid which plots near the olivine-plagioclase boundary originated by equilibrium 
processes governed by the phase relations of Fig. 1, i.e., by partial melting 
or crystallization, at pressures somewhat less than 10 kilobars, of a bulk 

(jJl'I'lifn 

system which lies - wi-t-h the (Mg,Fe).SiO.-SiO„-CaAl Si 0 system. If this model 
is assumed to be strictly valid then it is possible to "solve*' the problem with 
an approach described by Weill and McKay (1975) . Using distribution coefficients 
for Sr and the REE measured for plagioclase and olivine and assuming some initial value 
for the concentration of these elements in the parental material, the mass balance 



of plagioclase and olivine in the equilibrated solid assemblage. This solid 
assemblage, when added to the liquid C, constitutes the parental material. In 
the present calculations we have assumed an initial concentration of 5X chondritic 
levels for the LIL in the parental material, . Crystallization experiments in the 
pseudoternary system (Weill and McKay, 1975) indicate that liquid C is 
equilibrated with plagioclase and olivine at approximately 125 0°C , and we 
interpolated the plagioclase distribution coefficients of Table 2 to this temperature. 
The temperature dependence of olivine distribution coefficients is not known and 

was assumed to be similar to ^Q^/liq 

under strongly reducing conditions. 

The best fit (using least squares criteria applied to solutions of equation [1] 
for each of the trace elements) to these data yields a parental composition shown 
as C* in Fig. 5. If this material originally had 5x chondritic abundances of LIL, then 
approximately 4 wt % of liquid equilibrated with a solid assemblage dominated by 
olivine and plagioclase would have an LIL abundance pattern' closely approximating 
that of KREEP-rich composition C. There are many factors which prevent a rigorous 
assessment of such results. Among others, we can list the more obvious ones: 

(1) Modest departures of bulk composition as well as increased pressure can result in 


we .‘.have adopted the values at 1200°C. 





- 11 - 


appreciably different "phase diagrams" such that liquid C is not necessarily 
saturated with olivine and plagioclase as we have assumed. (2) The assumption that the 
LIL elements in the parental material have a primitive lunar (i.e.j flat 5X chondritic) 
abundance pattern is somev/hat arbitrary. (3) If the source of KREEP was itself 
formed as a direct or indirect consequence of fractionation in a '' lunar 

magma ocean, then there is always the possibility that the source mineralogy 
and bulk composition was dominated by cumulates while the trace element abundances 

Vte.ui> 

were influenced by varying amounts of intercumulus liquid. In-5^3few-of these and 
other uncertainties we believe that these calculations simply show that if the 
KElEEP-rich liquid precursors (e.g., liquid C) to the highly evolved igneous KREEP 
samples (e.g., 14163 and 15386) were derived from primitive lunar systems under 
conditions approaching equilibrium partitioning of trace elements, then that 
source material must have been rich in olivine and contained much smaller amounts 
of plagioclase. This is essentially the same general theme that emerged from our 
previous modeling (Weill and McKay, 1975) which considered the generation of 
evolved KREEP (for instance, near-peritectic liquids) by one-step partial melting 
of assemblages resembling the whole moon in bulk composition. 
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Eb-Sv isotopes 

The plagioclase/liquld distribution coefficients for Sr and Rb differ by 

almost two orders of magnitude (Table 2) . Any fractionation step in the evolution 

of the KREEP suite which involves plagioclase and liquid will therefore result in 

a change of Rb/Sr. We now consider whether the models we have proposed are consistent 

with the Rb/Sr isotopic data. If a closed .system inherited~a radiogenic strontium 

87 86 

content denoted by ( Sr/ Sr) at some time, T , before the present, then the 

o o 

present (^^Sr/^^Sr) and (^^Rb/^^Sr) are given by: 

P P 

(®^Sr/®^Sr) - (®^Sr/^^Sr) = (®\b/®Sr) (e^^o - 1) • [2] 

pop 

where X = 1.39 x 10 yr ^ is the ®^Rb decay constant. Only if a system has 

actually remained closed to Rb/Sr fractionation between T^ and the present 

does T^ really correspond to a time in the past when the system inherited 
87 86 

( Sr/ Sr)^, However, according to equation [2], if a system is assigned a 
87 86 

value of ( Sr/ Sr)^ , a value of T^ can always be calculated regardless of 


whether the system has actually undergone Rb/Sr fractionation or not. For a 
system which has suffered Rb/Sr fractionation, T^ has no real significance 
except to denote the time in the past when 



evolved without 


the system would have inherited (^^Sr/^^Sr)^ 'had 'it 

Rb/Sr fractionation. For this reason, Tq' is usually referred to as a model 

"age" and given the notation Tjj or Tjr where the subscripts signify model 

'87 86 

and assigned Initial ratio of ( Sr/ Sr) respectively. 

Schoenfeld and Meyer (1972) have shown that if a system actually inherited 
87 86 

(Sr/ Sr)^ at and subsequently underwent a Rb/Sr fractionation step at 
Tx , then the model age defined by equation [2] is given by 


Rl 


[3] 


where 


(Rb/Sr) 

(Rb/Sr) 

' parent 


Schoenfeld (1976) extended the same approach to consider’ two fractionation steps, 
and it can easily be shown that for three fractionation steps 



(e^^o - e^^l) (e^*^! - 


e^^2) ^ (e^‘^2 - e^'^3) 


+ e 


XT3 


E4] 




R2R2R3 



•14 


We now apply the formalism - o - tlino above to the models we have discussed for 

the evolution of the vatious members of the KREEP suite. In the calculations 

we have adopted values of = 4.60 b.y. and (^^Sr/^^Sr)^ = 0.69903 (LUWI of 

Nyquist et al. , )♦ Most of the whole-rock isotopic data we use were 

determined by Nyquist and coworkers. We have also used data determined by 

87 86 

Papanastassiou and coworkers and have adjusted ( Sr/ Sr) for interlaboratory 
bias. None of our arguments are significantly affected by changes of ±0.5 b.y. 

07 or 

in Tq or by adopting ( Sr/ Sr)^ = 0.69898 (BABI of Papanastassiou et al . , 1969). 

Composition C (Fig. 5) is the least evolved of the KREEP suite we have 

discussed and may be a precursor type for the other members of the suite. In 
I ■ . at/e»-ayes of F<^(>lished. 

Table j we have listed/|^L Rb/Sr isotopic data avm M gpn - for rocks that 

met the criteria for parental material of 14163 or the calculated P 

composition. The details of petrogenesis for liquids with compositions close 

to C are not known, but for the present discussion we need simply to adopt the 

view that compositions of this type were produced in a single differentiation 

step from primitive lunar materials. Equation [3] may be solved for Tjj^ vs. Rj^. 

The solution is shown in Fig. 6 . If we adopt the value of (Rb/Sr) = 0.0067 

parent 

as proposed for the whole moon by Taylor and Jakes (1974) , the time of the 

va\oe 

differentiation step is 4.44 b.y. This/\ls essentially equal to the 4.42 b.y. 

proposed by Tera and Wasserburg (1974) on the basis of U— Pb isotopes for the 

time of large scale differentiation of the moon. 

We have seen that liquids of cotectic composition C may be typical of the 
KfiEEP-rich wore 

precursors of/]highly evolved KREEP samples. such as 15386. According to 

our model,' 15386 may be the result of three differentiation steps; 

Ti T 2 T3 

primitive parent >- C — >■ P > 15386. The internal isochron age of 

15386 is 3.94 b.y. (Nyquist et al. , 1975) and the simplest interpretation is to 

07 Q£ 

equate this with T 3 , the time of the most recent differentiation. ( ^Rb/ °Sr) 

OBIGINAL PAGE IS 
OP POOR QUALITY ' 



values for the whole moon, compositions C and P and 15386 are all listed in 
1 

Table X» snd these values fix R 2 > ^ 3 * Furthermore, = 4.44 b.y. and 

T^ = 3.94 b-y., and we can solve equation [4] for the model age of 15386 vs. T 2 . 

The solution is shown in Fig. 7- The model* age of 15386 is 4.26 b.y., 

corresponding to T 2 = 4.41 b.y. This result is not at all sensitive to variations 

in ( Sr/ Sr)o. It is somewhat sensitive to variations in the last fractionation 

step, but it is not possible to place meaningful accuracy limits on R^. However, 

ore 

is calculated from data which -±a independent of the isotopic data, and we 

believe that the result T 2 - Tj^ is not fortuitous. The production of liquid C 

and the subsequent differentiation to liquid P in a short time is appealing from 

fn 

a physical point of view. It is difficult to store a liquid for long periods :©a the 

moon, and the alternative would require a later heat source to remelt C. 

The present-day isotopic ratios of P, C and 15386 are shown in Fig, 8 along 

those of some i^REEP-'r'icU 

with^ Apollo 14 and 15 crystalline KREEP samples and^the dSEERy rocks found to 
be possible parents for 14163 and composition P. Also shown is the isotopic 
composition of the whole moon (M) calculated from the whole moon Rb/Sr of 
Taylor and Jakes (1974), the (^^Sr/^^Sr)^ = LUNI of Nyquist et (1974) and 
Tq = 4.60 b.y. M, C, and Pj as well as the possible parents of Pj lie very close to 
the -4.44 b.y. isochron, while the Apollo 14 and 15 crystalline KREEP rocks are part 
of the 3,94 b.y. isochron defined by P and 15386. 

Another perspective of the Rb/Sr isotope evolution of these samples is shown 
in Fig. 9A where we plot the isotopic compositions of the samples 3.94 b-y. ago 
calculated on our model assumption that no Rb/Sr fractionation has occurred since 

that time. Most of the Apollo 14 and 15 crystalline KREEP samples have the same 

87 86 ^ 

C'Srr'^Sr) as P and 15386 within analytical error limits, consistent with the 

production of these evolved KREEP samples 3.94 b.y. ago by an igneous process 

87 86 

which involved variable enrichments of ( Rb/ Sr). The logical candidate for such 

0KIGMALPA6B® 

OF POOR QUALlTi 
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a process is partial melting of crystallized P or fractional crystallization of 

liquid P. The samples that were found to have compositions that cou/d parental 

to peritectic liquids are still on the isochron defined by M, C and P, indicating that 

these samples were not appreciably fractionated 3.94 b.y. ago. 

QS they hove besh 

In Fig. 9B we have plotted the isotopic compositions/\4. 44 b.y ago. The 

crystalline Apollo 14 and 15 KREEP samples are not included because they did not 

exist as such according to our working model. The precursors to KREEP 

87 86 

peritectic liquids all have ( ^Sr/ °Sr) identical to the whole moon and P. The 


, 87 86 * 

spread of ( ^Rb/ Sr) at this time is consistent with derivation of a suite of 
liquids from a primitive lunar composition. Note that it would be possible to 
generate a range of Rb/Sr values in such liquids by varying the relative propor- 
tions of liquid to crystalline material (most importantly, plagioclase) . Compo- 
sition P is just one possible member of the suite generated at this time. In general, 

compositions close to P would be equilibrated with more residual plagioclase 

87 86 

(and hence undergo greater increases of (°^Rb/ Sr) relative to M) than liquids 

further up/\the olivine-plagioclase cotectic. Also, since it is possible for a 

liquid of peritectic composition to coexist with variable amounts of plagioclase 

we could even expect that peritectic liquids (as defined by their major element 

composition) derived at this time might hav^ariable ( Rb/ Sr) . The composition 

P shown in our plot is just one of the series singled out because of its 

calculated parental relationship to 15386. If another member of the series 

generated 4.44 b.y. ago subsequently differentiated at 3.94 b.y., its daughters 

87 86 

would have inherited somewhat higher or lower ( Sr/ Sr) at that time (cf. the 

c. 

Apollo 14 and 15 crystalline KEIEEP in Fig. 9A) . 



SUMMARY 


K^EEP- rich 

The model which we have presented to account for a range of lunar 

rocks is summarized in Fig. 10. It relates much of the KREEP suite in a scheme 
of igneous differentiation which is consistent with major and trace element 
concentrations, the general features of experimentally determined phase 
equilibria and Rb-Sr evolution. The "proto-KREEP" from which more evolved 

L 

KREEP is derived strongly resembles a common highland rock type, the ^ow-K 
Fra Mauro basalts. Although the model describes a possible three-step 
evolution for the more evolved crystalline KREEP samples such as 15386, the 
most probable time scale according to the isotopic data only allows resolution 
of two igneous activity pulses. The formation of the highly evolved 15386 
occurred during igneous activity about 3.9 b.y. ago, at a time when the rate 
6f Impacts on the lunar surface had slowed down and mare basalt generation was 
beginning . 

The origin of the KREEP* precursors (LKFM(?) peritectic liquid types) 

at 4.4 b.y. is not well delineated but most probably involved liquid differentiation 

from a relatively primitive lunar source which contained much olivine and less 
' tbs po55ibiflty tifft 

plagioclase. We cannot rule out/] a limited amount of other refractory mineral phases 

of HR££P-rich 

oAso left behind after this liquid type was generated ./j the ; _ systems derived at 

this time do not seem to have been subsequently strongly disturbed in major, trace 
or' isotopic composition. 
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Table 1. Major element, trace element and Rb-Sr isotopic compositions of 


several members of the KREEP suite. The average highlands and 
whole Moon compositions proposed by Taylor and Jakes (1974) are 
given for comparison. See Figures 3 and 8 for sources of data. 


k 

15386 

14163 

pi 

c' 

LKFM^ 

HL^ 


wt% 








SiO^ 

50.8 

48.0 

49.1 

45.8 

46.6 

44.9 

44.0 

Ti02 

2.2 

1.8 

2.3 

1.3 

1.3 

0.6 

0.3 

A1,0^ 

14.8 

17.6 

15.9 

19.0 

18.8 

24.6 

8,2 

FeO 

10.6 

10.4 

8.2 

8.6 

9.7 

6.6 

10.5 

MgO 

8.2 

9.2 

11.8 

12.7 

11.0 

8.6 

31.0 

CaO 

9.7 

11.2 

12.1 

11.3 

11.6 

14.2 

6.0 

ppm 








Rb 

18.5 

13 

11 

5 . 6 

2.7 

1.7 

0.3 

Sr 

187 

180 

180 

167 

140 

200 

43 

Ce 

211 

200 

126 

78 

81 

22, 

3.1 

Sm 

37 

32 

22 

14 

13 

3.3 

0.8 

Eu 

2.7 

2.8 

2.2 

1.8 

1.8 

1.7 

0.3 

Yb 

24 

24 

15 

9.7 

9.8 

2.3 

0.7 

®’sr/«sr 

0.71640 


0.71049 

7 

0.7061^ 




®\b/«^Sr 

0.285 


0.180 

0.112 




"^LUNI 

4.26 


4.45 

4.47 






1. Peritectic liquid which may have been parental to 15386. 
Composition calculated from 15386 as explained in text. 
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Table 1. (continued) 

2. Average of samples found to be potential parents to P and/or 14163. 

3. Apollo 15 Low-K Fra Mauro basalt of Reid et al. (1972) . Trace element composition 
from Taylor (1973) . 

4. Proposed average composition of highland crust (Taylor and Jakes, 1974). 

5. Proposed composition of whole moon (Taylor and Jakes, 1974). 

6. = 0.69903 (Nyquist aT . , 1974). 
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Table 2. Experimentally determined crystal/liquid distribution coefficients 

(wt. basis) for ilmenite, plagioclase, Ca-poor pyroxene and olivine. 
Values in italics are from Weill and McKay(1975). Other values are 
from the present study. The uncertainty generally represents the 
standard deviation of the mean of a number of determinations. 
Representative compositions of coexisting phases are given in notes a-f. 


Ilm^ Plag^ Plag^ Px*^ Px® - 01^ 

(1140'C) (1340°C) (1200^0 (1340'^C) (1200'’C) (1200'’C) 



Oxide values are in weight %. To facilitate location in the pseudoternary 
diagram (Fig. 1), mole % silica (SI), anorthite (AN), and olivine (OL) are 
given for glasses other than mare basalt, 
a. Ilmenite; Ti 02 = 55 . 9 , A1^0^==0.37, FeO=34.5, Mg0=7.07, Cr 202 = 1.79 

Glass: Si02=42.4, Ti 02 = 10 . 2 , Al 202 = 7 . 99 , Fe0=19.0, Mg0=6.26, Ca0=8. 98 ; Cr^ 03 = O./O 
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Table 2. (continued) 


b. Plagioclase; Si02=44.7, Al20^=36.2, Fe0=0.29, Mg0=0.42, CaO=19.8 

Glass: SiO =48.0, TiO =.62, Al^O =23.1, Fe0=3.50, Mg0=10.1, CaO=11.6 (SI. ,AN„_OL„,) 
L 2. 2 6 Jd 6o 2b 


c. Plagioclase; Si02=45.3, Al20^=34.0, FeO=0.33, Mg0=0.44, CaO=18.7 

Glass; SiO =50.2, TiO =1.49, A1^0.=17.7, FeO=8.65, Mg0=9.12, CaO=9.27 (SI . „AN„ ,OL. .) 

z 2 2 6 4o zD zu 

d. Pyroxene: Si02=59.3, Ti02=0.05, Al202= 0.53, Fe0=4.34, MgO=36. 9,CaO=0. 26 

Glass: 'SiO-=57.4, Ti0.=.-60, Al.,0.,=10. 64, Fe0=6.92, Mg0=15.7, Ca0=5.-67 (SI,.„AN, ^0L„„) 

2 2 2 6 by 12 2y 


e. Pyroxene: Si0^=54.2, TiQ2= 0.3, Al202=5.3, Fe0=8.7, Mg0=30.4, Ca0=1.3 

Glass: Si02=48.9, Ti02=1.10, Al203=16.6, Fa0=9.37, Mg0=8,90, Ca0=9,23 (si^gAN^^OL^^) 

f. Olivine; 8iO„=40.3, Fe0=13.6, Mg0=46.2, Ca0=0.20 

Glass: SiO„=48.6, T10„=.93, A1„0 =16.7, Fe0=9.20, MgO=9.75, CaO= 9.04 (SI, ^AN^^0L^„) 
z z 2 6 46 2b 2y 

g. Values for f = 10“^^'^at 1200“C and'l0“^^‘ at 1340° C. See Fig. 2 of this work 

^2 

and Fig. 5 of Weill and McKay (1975) for oxygen fugacity dependence. All other 


distribution coefficients are independent of f^ over the range of fugacities 

studied, log f = (6. 57-27215 /T) ± 2. 

* 
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FIGURE CAPTIONS 

KR££P-rfcL 

Figure 1. Representation' of bulk compositions of 45REEPy materials from Apollo 
14-17 within (Fe,Mg) 2 S 10 ^-Si 02 -CaAl 2 Si 20 g pseudoternary liquidus 
diagram. Average lunar highlands composition of Taylor and Jakes(1974) 
is shown for comparison. 


Figure 2, Variation of Cr distribution coefficients (weight basis) for 


olivine/liquid .and orthopyroxene ./liquid with .oxygen fugacity. Error 

bars represent twice the standard deviation of the mean of a number of 

replicate determinations at each oxygen fugacity. 

Figure 3. Calculated Eu contents and Mg/ (Mg+Fe) of hypothetical peritectic liquids 

e 

having the same Sm content as 14163 . Each point represents the liquid 



produced by differentiation of a different highland sample. See text 

for details of the calculation. Actual Eu content and Mg/ (Mg+Fe) 

of 14163 are indicated. Estimated uncertainties (±30% for Eu, ±20% for 

Mg/ (Mg+Fe)) in calculated values for hypothetical liquids are indicated 

poten /-/a-£ 

by the rectangle around 14163. Highland samples which are poDciblc - 
parents to 14163 must produce liquids which plot wi.thin this rectangle. 
Liquid produced from the average highlands composition of Taylor and 
Jakes(1974) is indicated. Sources for data: Drake et al.(1973), Haskin 
^^.(1973), Hubbard ^ (1972,1973), Laul ^ al. (1972) , LSPET(1972, 

1973) ,Nava(1975) , Nakamura ^ al . (1973) , Philpotts ^ al . (1972,1974) , 
Rhodes and Hubbard (1973) , Rhodes et al. (1974) , Rose et al.(1972), 

Taylor et al . (1972,1973) , Wanke et al. (1972) , and Wiesmann and Hubbard 
(1973) . 


Figure 4. (a) Pseudoternary projection of bulk compositions of highland samples 
<tS pofenfia^. po.ireiJs 

tested for — pe -»& iblg r paternity - to 14163 (b) Bulk compositions 

- poi&n ha£ 

of s 0 japles found to be po ss ib l e parents to 14163, Le.^ samples generating 
liquids which plot within the rectangle surrounding 14163 in Figure 3. 



Figure 5. Pseudoternary representations of important compositions in our model for 


the evolution -of KREEP. P is the peritectic composition which differentiated 

• Campos/h^’*- 

to form 15386. 14163 is another example of a KEEEP-rich peritectic - liquid -. 

po*4e>^ 

C is the average composition of all highlands samples which are - poss - ih -L&- 

posS‘l>te. 

parents to P or 14163. C' is the composition of -the/] parent material to C, 
Ca.lcufa.'f'cd e.7: p/oj ned m {-he. 


— te - C had 5x chond r ltic LIL -abundanc es-. WM -is the whole .moon composition 

V 

proposed by Taylor and Jakes (1974). 

Figure 6. Solution of Equation [3] for T^, the time of formation of composition 
C (assuming a one-step process)jVS. the amount of Rb/Sr fractionation 


Figure 7. 


Figure 8. 


I'l 

P 


occurring during its formation. The dotted vertical line is the Rb/Sr 

dS 

fractionation if C was formed from material with Rb/Sr>^ proposed by 
Taylor and Jakes (1974) for the bulk moon. The dotted horizontal line is 
the resulting value of T- . 

Solution of Equation [4] for T^, the time of - f orm a t -i-o n - of composition P 
from^C, vs. the model age of 15386. If the bulk moon Rb/Sr of Taylor and 
Jakes(1974) is adopted, R^, R 2 , and R^ are fixed by the data in Table r. 

//t- 

T^ is determined fey Figure 6, and T^ is assumed to be 3.94 AE, the 
internal isochron age of 15386 (Nyquist ^ a^. , 1975) . The dotted 
vertical line indicates the model age of 15386 reported by Nyquist 
et al . (1974) , and the dotted horizontal line represents the resulting 
value of T 2 * T^ is indicated for comparison. 

Rb-Sr evolution diagram showing presently observed or calculated 

isotopic ratios for 15386, P, C, and the whole moon (M) . According to our 

model, M, C, and P define a 4.44AE isochron, while P and 15386 define 

Q^so 

a 3.94AE isochron. Isotopic ratios are^shown for a number of highland 
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samples which are ^oos- i - h - le parents to P (small p's) and a number of 
Apollo 14 and 15 crystalline KREEP samples (x's). Sources of data: 
Nyquist et al,>(1973 ,1974,1975) , Papanastassiou and Wasserburg(1971) , 
and Wasserburg and Papanastassiou (1971) . 

Figure 9. Rb-Sr evolution diagrams showing isotopic ratios at time T, assuming 
closed system evolution from T until the present. Symbols are the 

C.rys.-faf/rlLa.-^rts<^ 

same as in Figure 9. (a) T=3.94AE, the prof - er for the- f orm atlox 

of 15386 (b) T=4.44AE, the time we calculate for the formation 
of C and P from M. 

Figure 10, Schematic representation of the model we propose for the origin of 


KREEP and the relationship between different materials. F is 
the fraction of liquid resulting from <any differentiation step relative 
to the immediate parent material for that step. 
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The partitioning of Mg, Fe, Sr, Ce, Sm, Eu, and Yb 
in lunar igneous systems and a possible origin of KREEP 
by equilibrium partial melting 

D, F. Weill, and G A McKay 

Department of Geology, University of Oregon, Eugene, Oregon 97403 

Abstract — ^The solid/liquid distnbution of Mg, Fe, Sr, Ce, Sm, Eu, and Yb have been determined over 
erystallization paths in the system (Mg, Fe)3Si04-Si02-CaAl2Si208 Varying degrees (1-5 wt%) of 
partial melting of an ohvine + orthopyroxene + plagioclase rock yield liquids with large-iron hthophile 
(LIL) and major-element concentrations strikingly similar to those of some Apollo 15 and 16 KREEP 
rocks The parent assemblages from which these KREEP melts can be derived have a bulk composition 
s imil ar to thatpreviously proposedfor the whole moon Partial melts derived from a sourcewith the same 
mineralogy but with a chemical composition corresponding to the average lunar highland do not 
simultaneously satisfy the major-element and LIL concentration characteristics of KREEP 


Introduction 

The differences of chemical composition observed m lunar rocks can either 
be attributed to heterogeneities associated with the original accretion of the moon 
or to subsequent large scale solid-liquid phase separation processes Geochemical 
arguments (Brett, 1973, Duncan et aU 1973; Taylor and Jakes, 1974) indicate a 
relatively homogeneous accretion. Furthermore, rapid accretion (Mizutani, 1972) 
could have resulted m an early liquid outer shell extending to a depth of at least 
several hundred kilometers. Cooling and crystallization of this initial liquid shell 
would result m chemical diiferentiation that would have characteristics imposed 
by solid-liquid partitioning of the elements. Numerous models have been prop- 
osed to explain how an initial global igneous complex on the moon might account 
directly (or indirectly through subsequent internal partial melting) for the main 
rock types observed m the lunar sample collection Volcamsm and large impacts 
on the lunar surface could have served to bring representative rock types or their 
partially molten derivatives to the surface It has also been proposed (Warner et 
at, 1974) that surface impacts may themselves cause differentiation of the lunar 
crustal material via liquid-solid separation after impact melting. 

In one form or another, almost all models for the geochemical evolution of the 
moon relate the main lunar rock types to each other via a relatively simple 
sequence of crystallization and/or partial melting processes One of the principal 
tests of such models is a verification of compositional trends by means of 
experimental crystallization studies Since both major- and trace-element con- 
centrations in the lunar samples place critical constraints on the models, it is 
important that experimental verifications include data on both types of con- 
stituents. We have previously described some of the important factors to be 
considered m determining the evolution of elemental distributions during lunar 
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Igneous differentiation processes (Weill et al, 1974) and applied them to the 
evolution of Sm (indicative of “typical” REE) and Eu (“anomalous” REE) 
concentrations dunng various postulated fractional crystallization and partial 
melting sequences in lunar igneous systems. 

The solid/liquid distribution coefficients we employed for Sm and Eu (Weill 
and Drake, 1973; Sun et al, 1974, Grutzeck et al , 1974; Drake and Weill, 1975) 
were necessarily approximate because they had been determined experimentally 
m systems which departed in composition from the lunar systems to which we 
applied them Furthermore, our postulated crystallization sequences were based 
entirely on the olivine-plagioclase-silica pseudoternary hquidus diagram (Walker 
et al , 1972, 1973) which does not contain any information about the interphase 
distribution of the major elements Mg and Fe Nevertheless, within the limits 
imposed by the data it was shown that the evolution of liquids with KREEP 
characteristics could most easily be accomplished by partial melting and that there 
was a rather tight mutual constraint between the assumed Sm and Eu abundances 
of the parent rock, the mineralogy of the refractory parent assemblage, and the 
degree of partial melting 

We have now added to our experimental data base by (a) Directly determining 
the crystallization sequences of a molten average highland composition (Taylor 
and Jakes, 1974) and a much more ferromagnesian liquid These experiments yield 
information on the distribution of the major elements. Mg and Fe, between 
olivine, orthopyroxene and coexisting liquids, and thus allow us to combine the 
constraints imposed by the evolution of these important major-element indicators 
with those imposed by the trace elements, (b) Directly measuring the solid/liquid 
distribution coefficients of the trace elements Sr, Ce, Sm, Eu, and Yb for 
plagioclase and orthopyroxene over the same crystallization sequences. This 
extends our data to five trace elements, including a control on the relative 
enrichment of light versus heavy REE Also, since the distribution coefficients 
have now been measured directly m the system to which the modeling will be 
applied, the compositional dependence is no longer a source of uncertainty. 

fn this paper we report the new experimental data, and we present an analysis 
of the possibility of generating KREEP-type liquids by one step partial melting or 
crystallization of systems with whole-moon abundances of LIL trace elements 
and a subsolidus mineralogy dominated by olivine, orthopyroxene, and plagio- 
clase. We also discuss the plausibility of deriving KREEP-type liquids by impact 
melting differentiation of the average highland surface. 


Experimental Methods 

The average highland (HL) and ferromagnesian. (FMG) crystallization sequences were experimen- 
tally determined by holding the appropriate bulk compositions at temperature and oxygen fugacity (Hi 
and C02 gas mixtures) in SiC resistance furnaces for penods of 24-^8 hr Charges were held m Pt-Rh 
wire loops formed from 2 5 x 10"^-m diameter wire which had been previously held m contact with the 
same molten composition at 1400”C and the I-W equilibrium oxygen fugacity for 48 hr The runs 
themselves were ail carried out at oxygen fugacities for which the Fe’TFe** m the liquid is negligible 
Compositions of the individual phases are based on a minimum of 10 point analyses per phase The 



The partitioning of Mg, Fe, Sr, Ce, Sm, Eu, and Yb m lunar igneous systems 


1145 


modes were determined from the bulk composition and the individual phase compositions by 
multilinear regression analysis After the liquid and solid compositions at various stages of the 
crystallization sequences were determined as described above, new experimental charges were 
prepared with a bulk composition equal to the liquid plus a small amount of plagioclase and/or 
pyroxene These compositions were doped with an additional 1-3 wt % of one of the LIL “trace” 
elements and held at temperature and oxygen fugacity for 100-150 hr After quenching, solid and glass 
were analyzed for the LIL element (distnbution coefficient) and major constituents (to check that the 
major-element compositions of all phases were close to those determined previously for the 
crystallization sequence) In effect, the trace-element distribution coefficients are determined in bulk 
systems that he slightly off the multiple saturation boundaries as depicted, say, m the ol-opx-plag 
pseudoternary liquidus system Each LIL and each sohd-liquid pair is studied in a separate group of 
experiments in order to circumvent problems arising from excessive enrichment in “trace” elements 
and from working with liquid-poor, multiply-saturated systems Runs used- to determine LIL trace 
element distribution coefficients were performed over a range of oxygen fugacities from “air” to two 
orders of magnitude below the I-W equilibrium curve Only the Eu plagioclase/liquid distribution 
coefficient exhibits significant dependence on oxygen fugacity 


Experimental Results 

The results of the experimental crystallization runs are given in Table 1 The 
Fe and Mg solid/liquid distribution coefficients for olivine and orthopyroxene are 
plotted against reciprocal temperature in Fig 1, and the linear regression 
coefficients are given in Table 2. Crystallization sequences for new bulk composi- 
tions can be calculated from the experimentally determined sequences, effectively 
stretching the applicability of the experimental set. For instance, the addition (or 
subtraction) of CaAUSUOs to bulk composition HL simply results in an equivalent 
increment of plagioclase at each step of the new crystallization sequence The 
addition of (Mg, Fe) 2 Si 04 to bulk composition FMG similarly results in additional 
olivme, but because of the Mg-Fe exchange between olivine, liquid, and or- 
thopyroxene during equilibrium crystallization the new crystallization sequence 
can be approximately calculated only with a knowledge of the Mg and Fe 
distnbution coefficients. In Table 3 we list a crystallization sequence calculated 
for a liquid with a bulk composition corresponding approximately to the liquid at 
1380°C in the original FMG sequence (Table 1) plus 54 wt % of the olivine (F 093 ) 
with which this liquid is in equilibrium This calculated sequence is referred to as 
“FMG ” Various aspects of the crystallization sequences HL and “FMG” are 
shown in Figs 2, 3, and 4 

The distribution of Sr, Ce, Eu, Sm, and Yb between plagioclase, orthopyrox- 
ene, and liquid is shown in Table 4 These coefficients are valid for the solid— liquid 
pairs which approximate in composition the coexisting phases along the crystalli- 
zation paths, i.e , near the boundary lines and reaction point of the 
olivme-plagioclase-silica pseudoternary system To the extent that their composi- 
tion dependence is not pronounced within this ternary system they are roughly 
applicable to coexisting phases within it Notice that although the coefficients are 
given at two temperatures, both the temperature and composition dependence are 
inseparably incorporated in the data. Low-pressure melting of an olivme + 
plagioclase + orthopyroxene assemblage begins slightly below 1200°C, and the 
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Table 1 Experimental crystallization of average highland (HL) and ferromagnesian (FMG) composi- 
tion liquids 


T (TO 

Phase (wt %) 

SiOg 

TiOi 

AlzOa 

CraOs 

FeO 

MgO 

CaO 

1400 

HL hq(lOO) 

45 2 

0 63 

25 3 

010 

65 

80 

14 2 

1320 

liq (77) 

45 7 

0 78 

22 2 

0 14 

S3 

103 

12 6 


plag(23) 

441 

0 00 

360 

0 00 

0 28 

0 37 

193 

1280 

liq(67) 

45 8 

' 0 87 

20 4 

0 11 

9.4 

117 

11 8 


plag(33) 

43 6 

0 00 

35 9 

005 

0 45 

0.41 

19 6 

1260 

hq(58) 

46 4 

1 10 

18 9 

0 15 

10 5 

118 

112 


plag(40) 

43 5 

0 02 

36 0 

0 01 

045 

0 37 

197 


ol(2) 

40 5 

0 00 

026 

0 10 

12 2 

46.4 

0 58 

1220 

hq(38) 

48 0 

1 64 

16 8 

0 18 

12 8 

10 0 

10 6 


plag(53) 

44 2 

0 02 

35 5 

0 02 

0 43 

041 

194 


ol(9) 

39.6 

0 02 

022 

0.16 

15 9 

43 6 

0 52 

1200"= 

hq-l-plag-4-opx-fol 








1400 

FMG hq(lOO) 

51 1 

031 

11 85 

022 

923 

21 1 

61 

1380 

hq(93) 

514 

0 33 

12 9 

021 

— 

190 

6.6 


0l(7) 

42 1 

0 00 

0 01 

0.12 

59 

516 

0 24 

1320 

hq(79) 

52 3 

0 39 

14 8 

0 22 

97 

15.1 

75 


oI(13) 

412 

0 00 

0 07 

017 

87 

49 6 

0-26 


opx(8) 

57.1 

0 00 

1 20 

0 32 

53 

35 7 

044 

1280 

hq(72) 

52.0 

0 42 

161 

019 

97 

13 4 

83 


0l(13) 

410 

0 00 

0.09 

0.15 

101 

48 4 

027 


opx(15) 

56 1 

0 03 

1.93 

0 45 

65 

34 2 

0.78 

1260 

Iiq(61) 

50 8 

0 48 

17 7 

013 

9.9 

116 

93 


ol(lO) 

406 

0 00 

0 06 

0.11 

11 3 

47 7 

024 


opx(29) 

54 9 

0 06 

3 7 

0 50 

7.2 

32 6 

1.07 

1240 

hq(58) 

50 2 

0 50 

17.9 

010 

10 3 

113 

98 


ol(6) 

40 5 

0.00 

— 

0 09 

11 2 

48 0 

0 21 


opx(36) 

55 0 

0 08 

40 

0 48 

71 

32 3 

1 11 

1220 

hq(54) 

50 4 

0 49 

182 

0 09 

10 0 

10 4 

10.4 


oI(6) 

40 7 

0 00 

006 

008 

12 4 

46 5 

0 28 


opx(40) 

54 1 

0 12 

50 

0 42 

77 

31.3 

129 

1200’* 

liq + ol + plag + opx 









*Low-quality analyses sufficient only for phase identification 


Table 2 Solid/liquid (weight basis) distnbution coefficients for Mg 
and Fe in olivine and orthopyroxene along the crystallization paths 
FMG and HL Coefficients in the equation, logiol? — A +B/T, 
where T is in °K The Fe-Mg exchange coefficient as defined by 
Roeder and Emslie (1970), = (X^,,)(Xi,J/(X;.c)(XM^, is equiv- 

alenttoDpc/J^Mg 


Solid 

Element 

A 

B 

Olivine 

Mg 

-1 451 

3138 

Olivine * 

Fe 



Orthopyroxene 

Mg 

-1232 

2557 

Orthopyroxene 

Fe 

-2 269 

3225 
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1220 1240 1260 1280 1300 1320 °C 



Fig 1 Fe and Mg distribution coefficients (weight basis) Open circles represent data 
from charges of bulk composition HL, closed circles from composition FMG Linear 
regression coefficients are given in Table 2 


melt does not vary greatly in composition provided the degree of partial melting is 
not large Consequently, in considering the partitioning of these trace elements 
between a partial melt and a refractory residue containing plagioclase and 
orthopyroxene, the coefficient values at 1200°C are representative It can be seen 
that plagioclase preferentially incorporates the light versus the heavy REE while 
orthopyroxene discnminates even more strongly in favor of the heavy REE 
Orthopyroxene displays very little potential for creating anomalous Eu 
abundances. In Fig 5 we have plotted the variation of the Eu distribution between 
plagioclase and liquid as a function of oxygen fugacity, showing the very large 
potential for anomalous Eu abundances developed in plagioclase under strongly 
reducing conditions. The coefficients corresponding to lunar oxygen fugacities 




Wh Fraction 
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Table 3 Calculated crystallization sequence “FMG.” 



"FMG” HL 



864202468 


F {Wt Fraction Residual Liquid) 

Fig 2 Phase assemblages (weight fractions) during eqmlibnum crystallization sequ- 
ences HL and “FMG ” Assemblages at F = 0 correspond to the normative proportions 
for the bulk compositions, with Fe and Mg partitioned between olivine and orthopyrox- 
ene according to the distnbution coefRcients in Fig 1 
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Fig 3 Equilibrium crystallization path m the pseudoternary system 
(Mg, Fe) 2 Si 04 -Si 0 j-CaAl 2 Si 20 a for bulk composition HL Multiple-saturation curves of 
Walker et al (1973) are given for reference For each data point the temperature, mg 
value of the liquid, and the Fo and En contents of any coexisting olivine and 
orthopyroxene are given Crystallization below 1220°C leads the liquid to the four-phase 
peritectic "point,” where crystallization proceeds (along with dissolution of olivine) to 
completion The exact location of the end point vanes with the mg value and the 
minor-element content of the bulk composition, but lies within the indicated envelope 
The mg, Fo, and En at the end of crystallization were calculated using the distnbution 

coefficients m Fig 1 


Table 4 Plagioclase/hquid and orthopyroxene/hquid distribution coefficients for Sr, Ce, Sm, Eu, and 

Yb 




1340X 

Plagioclase 

1200°C 



Orthopyroxene 

1340“C 

1200°C 

- 


D 

s„ 

n 

D 

Sm 

n 

D 

5™ 

n 

D 

Sm 

n 

Sr 

1 53 


6 

1 63 


4 



2 

0 018 

003 

3 

Ce 

IB 


5 



5 



3 

0 009 

002 

5 

Sm 

mB 








2 

0 022 

001 

9 

Eu 

1 1* 

— 

— 


— 

— 



3 

0 022 

002 

7 

Yb 



10 



6 



2 

0 170 

010 

6 


■^cf Fig 5 
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Fig 4 Calculated equilibnum crystallization path for bulk composition “FMG ” Sym- 
bols same as in Fig 3 


(WeilL et al , 1974) have been incorporated in Table 4. The differences between the 
present set of plagioclase/hquid distribution coefficients and those determined by 
Drake and Weill (1975) are illustrated in Fig 6 These differences are indicative of 
the variation in the coefficients to be expected as a result of compositional 
variations in the two host phases (lunar *‘olivine-plagioclase-silica” liquids versus 
terrestrial andesites and andesitic basalts, and Anioo versus Ansj-Aiiss plagioclase). 

Discussion 

The crystallization of HL and FMG liquids are seen m Figs- 3 and 4 to conform 
quite closely to the predictions based on the olivine-plagioclase-sihca pseudoter- 
nary This may be interpreted as one more confirmation that the pseudoternary 
diagram is a realistic approximation in the sense that the primary crystallization 
“fields,” three-phase “boundaries,” and four-phase “points” are useful references 
for approximating phase proportions and compositions (with the exception of 
MgO and FeO) These exceptions turn out to be of interest in discussing the 
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Fig 5 Vanation of Deu for plagiDclase/liquid with oxygen fugacity 


possible origins of KREEP compositions by either partial melting or near- 
equilibrium crystallization. The major-element compositions of KREEP rocks 
typically plot in the area around the reaction “point” where the liquid is multiply 
saturated with olivine, orthopyroxene, and plagiodase This is compatible with an 
ongin of KREEP by relatively low degrees of partial melting as pointed out by 
Walker et ai (1972)'> KREEP has relatively “primitive”* mg values (mg == 




1152 


D F Weill and G A McKay 



Fig 6 Plagioclase/Iiquid and orthopyro\ene/Iiquid distribution coefficients Coefficients 
determined by Drake and Weill (1975) in terrestrial basaltic andesites are shown for 

companson 


100 Mg/(Mg + Fe), atomic basis). The relatively primitive mg and the highly 
differentiated LIL abundances are part of the puzzle of KREEP. The extensive 
fractional crystallization or very low degree of partial melting required to achieve 
the LIL enrichment have been considered hkely to produce mg values too low to 
be compatible with KREEP In Table 5 we list analytical data for three Apollo 16 
(mg = 69, 65, 66) and one Apollo 15 (mg = 58) KREEP -type rocks Referring to 
Figs. 2 and 3 we see that for the HL sequence this range of mg values is 
encountered only with a relatively large fraction of melt (> 30 wt %) such that the 
resulting enrichment {liquid/average highland) m the LIL elements would be far 
short of those typical of KREEP. This seems to be a problem m interpreting 
KREEP as the product of impact melting induced differentiation of the highland 
surface (Warner et ai, 1974) Referring to Figs- 2 and 4 we see that the partial 
melting of a more mafic material as typified by “FMG” would produce the correct 
major element concentrations (i e in the vicinity of the pseudo reaction point and 
with suitably high mg values) at low degrees of partial melting Some of the 
problems associated with producing KREEP by partial melting of plagioclase-rich 
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Table 5 Partial analyses of KREEP rocks Sources Hubbard ctal (1973), 
Rhodes and Hubbard (1973), Hubbard et al (1974), Wiesmann and 

Hubbard (1975) 


Wt% 

64815 

62235 

65015 

15386 

S 102 

47 0 

47 1 

47 2 

50 8 

AbO, 

17 6 

IS 9 

23 2 

14 8 

FeO 

95 

94 

87 

10 6 

MgO 

11 6 

10 0 

95 

82 

CaO 

11 9 

11 5 

101 

97 

ppm 

Sr 

144 

161 

164 

187 

Ce 

93 7 

153 

125 

211 

Sm 

17 6 

27 1 

22 2 

37 5 

Eu 

I 68 

2 03 

191 

2 72 

Yb 

12 6 

18 7 

15 3 

24 4 

mg 

69 

65 

66 

58 


rocks were previously discussed m terms of Sm and Eu abundances (Weill et aL, 
1974) Briefly, it was pointed out that a high concentration of plagioclase in the 
refractory residue (e g as m the highland surface) results in excessive Eu 
depletion relative to the overall enrichment m REE (e g Sm) found m KREEP 
The results discussed above tend to confirm the unsuitability of HE as a parent 
material m terms of some of the major-element requirements of KREEP. 

In the last analysis (no doomsday pun intended) all geochemical tests of 
models for the origin of lunar rock types such as KREEP reduce to plausibility 
arguments It must be admitted, however, that there is some arbitrariness about 
the selection of an ad hoc “FMG” composition which then “turns out” to be 
compatible with a KREEP partial melt Consequently, we must refrain from 
claiming too much for it, and simply be content to say that, if KREEP was formed 
by a one-step, near-eqmlibrium partial melting process, a parent material similar 
to “FMG” is not ruled out, while one such as HL is. 

With a knowledge of the distribution coefficients for five trace elements (and a 
few assumptions) it turns out that we can approach the problem in a more 
systematic way. If we assume that partial melting of a mineral assemblage 
dominantly composed of olivine, orthopyroxene, and plagioclase is the origin of 
KREEP, the distribution coefficients can be used to test whether this is possible 
and, if so, to calculate the refractory assemblage in terms of the proportions of the 
minerals Obviously, we must also assume that melting occurs at sufficiently 
modest pressures so that the experimental distribution coefficients, obtained at 
atmospheric pressure, are representative. Rough estimates of the differences in 
partial molar volumes (solid versus liquid) for the trace elements in question do 
not indicate any unusually strong pressure dependence. It would be nice to know 
the pressure dependence of the distribution coefficients, but, in the meantime, it 
seems reasonable to assume that they will not change greatly at a few kilobars. 
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The relative values (even more critical) should hold to even higher pressures 

The problem is best approached by considering the simple mass-balance for a 
trace-element partitioned between a partial melt and olivine, orthopyroxene, and 
plagioclase: 

D = Do, Wo, + Do,. Wop. + Dp, Wp, - (Co - CrF)/(l - F)Cl, 

where D is the composite sohd/liquid distribution coefficient, C„ is the bulk 
system concentration, Cl is the concentration in the melt (KREEP), F is the mass 
fraction of liquid, and Wo,, Wopx, and IV'pi are the mass fractions of the minerals 
normalized to a sum of unity. The models we present assume Co to be equal to 7 x 
chondntic abundances, and Cl is set equal to the concentration in a particular 
KREEP rock. We consider that partial melting occurs at a temperature of 
approximately 1190°C and Do,, and Dpi are taken from the data in Table 4 
extrapolated slightly from 1200°C to 1190'’C. Do, for all the LIL elements in 
question is assumed to be negligible (< 0.01) At any selected degree of partial 
melting (F) the above equation is linear in two independent variables (2 out of the 
3W’s). With five constraining trace elements we can obtain the best solution with a 
multilinear regression analysis at each selected value of F. A typical result of this 
procedure is shown in Fig. 7 where we plot the quality of the solution (r.m.s.d. 
between calculated and assumed D) against F A best-fit for rock 64815 is clearly 
indicated at F = 05. An element-by-element comparison of the calculated partial 
melt with 64815 is illustrated in Fig 9 on the familiar abundance pattern plot The 
mode of the refractory residue (i.e. Woi, Wapx, and Wpi) is plotted in Fig 8. The 
best solutions for other Apollo 16 KREEP-rich rocks, 62235 and 65015, are very 
similar, and, for these, we show only the calculated refractory modes in Fig 8. 
Apollo 15 rock 15386 is another, and perhaps better candidate for igneous 
KREEP, and the best solution for this rock (F = .01) is shown in Figs. 8 and 9. The 
refractory modes (approximately equivalent to the parent material for low degrees 
of partial melting) cluster rather tightly, indicating that these four KREEP partial 
melts could be produced by different degrees of partial melting of essentially the 
same parent rock type The proportions of plagioclase in the preferred parents 
(20-30 wt.%) IS slightly greater than that previously deduced solely on the basis of 
Sm and Eu abundances (Weill et al., 1974). The difference is in part due to the use 
of somewhat different distnbution coefficients As previously discussed', the 
present set is more appropriate 

The above approach can also be used assuming that the LIL concentrations m 
the KREEP source are equal to those of the average highlands (Taylor and Jakes, 
1974) As expected from our previous discussion of the HL composition, the 
resulting solutions yield parent materials that differ substantially from the average 
highlands in mg values. 

We have assumed LIL concentrations in the KREEP source that are not 
greatly different from recent estimates for the whole moon (Taylor and Jakes, 
1974), but we have assumed nothing about the relative proportions of minerals m 
the source Therefore, it is interesting, and not a product of circular reasoning, 
that this mineralogy results m a major-element composition that resembles some 
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F 

Fig 7 Quality of fit (root mean square of the percent difference between calculated and 
assumed D for five trace elements) of the calculated refractory mineral assemblages m 
equilibrium with a hquid of KREEP composition 64815 at varying degrees of melting 

Tabic 6 The bulk compositions of olivine + orthopyroxene +plagiocIase assemblages 
(7 X chondntic LIL abundances) which form partial melts corresponding to various KREEP 
samples For comparison, two estimates of the composition of the whole moon are also hsted. 



64815 
F= 05 

62235 
F= 02 

65015 
F= 03 

15386 
F= 01 

Whole 
moon (I) 

Whole 
moon (2) 

Sl02 

462 

45 1 

46 0 

43 9 

44 1 

42 6 

AhO, 

10 3 

90 

92 

70 

82 

82 

FeO 

75 

85 

82 

12 0 

10 5 

10 9 

MgO 

317 

32 4 

31 6 

33 2 

312 

30 6 

CaO(3) 

58 

50 

50 

39 

60 

76 

mg 

88 

87 

87 

83 

84 

83 


(1) Taylor and Jakes (1974) (2) Ganapathy and Anders (1974) All normalized to 100 wt % 
(3) In converting from the mineral assemblages shown in Fig- 8 to the bulk compositions listed 
here it was assumed that the pyroxene did not contain any CaO The low values of CaO 
. relative to whole Moon estimates are partly attnbutable to this procedure 




ppm /chord 


1156 


D F Weij^l and G A McKay 



Ol PI 


Fig 8 Refractory residual assemblages in equilibrium with four KREEP liquids LIL of 
bulk system IS assumed to be 7 x chondntic 



Fig 9 Comparison of calculated trace-element abundances (lines) with measured 
abundances (points) for KREEP samples 64815 and 15386 Liquids are in equdibnum 
with refractory residues whose modes are shown in Fig 8 Abundances were calculated 
as discussed m text, using the preferred F values indicated in Table 6 
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recent estimates for the whole moon In Table 6 we list the bulk composition of 
the parent material (refractory mineral residue plus the small fraction of KREEP 
melt) for each of the four preferred solutions. These are compared with two 
recent estimates of the whole-moon bulk composition. The agreement is not 
perfect, but, m view of the total range of compositions possible m the 
olivme-orthopyroxene-plagioclase system, we will adopt the classic stance and 
“stress the similarities rather than the differences ” It might also be stressed that if 
the whole-moon composition was ever a liquid as imphed in the primitive liquid 
shell model, then a mineral assemblage corresponding approximately to the 
undifferentiated composition might well be expected to occur near the surface 
We have presented experimental data on phase equilibria and the partitioning 
of two major elements (Mg and Fe) and five trace elements (Sr, Ce, Sm, Eu, and 
Yb) in lunar igneous systems and have attempted to show how such a combination 
of data provide a powerful constraint on geochemical modeling Obviously, each 
element and phase added to the list will provide much greater credibility to the 
decreasing number of models that remain plausible. There are many elements and 
phases left, and we can close on an upbeat note* more work can certainly be done 
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The effects of temperature, and composition on the electncal conductivity of silicate liquids have been ex- 
perimentally determined from 1200 to 1550°C under a range of /Qj conditions sufficient to change the oxidation 
state of Fe from predominantly Fe’"*" to Fe^"*" Oxidation of ferrous to ferric iron in the melt has no measurable ef- 
fect On the conductmty of melts with relatively low ratios of divalent to univalent cations. Under strongly oxidizing 
conditions a minor decrease of conductivity is detected in melts with high ratios. It is concluded that 

for purposes of estimating the conductivity of magmatic liquids, /Oj may be ignored to a first approximation Both 
univalent and divalent cation transport is Involved m electrical conduction Melts relying heavily on divalent cations 
for conduction, i e melts with relatively large SM*'*'/2M‘'' ratios, show strong departures from Arrhenius tempera- 
ture dependence with the apparent activation energies decreasmg steadily as the temperature increases Conductivi- 
ties dominated by the univalent cations, in melts with relatively small ratios, show classical Airhemus 

temperature dependence. These observations are discussed in terms of the general characteristics of the melt struc- 
ture. 

Compositional variations withm the magmatic range account for much less than an order of mapitude variation 
in electrical conductivity at a fixed temperature Thu observation, combined with previous measurements of the con- 
ductivity of olivine (A. Duba, H.C Heard and R Schock, 1974) make it possible to state with reasonable confidence 
that melts occurring within the mantle will be more conductive by 3-4 oidets of magnitude than their refractory 
residues. Potential applications to geothermometry are discussed 


1. Introduction 

Experimental measurements of electrical conduc- 
tion in silicate liquids are of interest to the geoscien- 
tist for several reasons. High electrical conductivity 
anomalies in the mantle associated with a low seismic 
velocity zone are not compatible with the much lower 
conductivities of olivine recently measured by Duba 
et al. [1], and are likely to be related to the presence 

* Present address; Department of Geology and Geophysics, 
Yale University, New Haven, Connecticut 06520. 


of a more conductive silicate liquid phase. Waff [2] 
considered the problem of modeling the bulk electri- 
cal conductivity of a solid + liquid mantle. A general 
result of this study was that if there is nearly total 
liquid bndgmg along gram boundanes or edges and if 
the conductivity of the liquid is much greater than 
that of the solid phases, the bulk conductivity is deter- 
mined by the concentration of the hquid and its con- 
ductivity. Since the composition, temperature pres- 
sure and oxidation state of partial melts can be ex- 
pected to vary considerably, the effect of these vari- 
ables must'be established before it can be assumed 
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that > Osoiid under all probable mantle condi- 
tions. Such measurements also provide a direct means 
of studying the kinetics ol ionic transport.and inter- 
preting them in terms of melt structure under con- 
trolled (T, X, fQ^ conditions Conductivity measure- 
ments can be made with high precision and are there- 
fore potentially sensitive to physical and chemical 
changes in the melt 

Measurements of electrical conductivity of molten 
and partially molten rock systems were first reported 
by Barus and Iddings in 1882 [3] and have been con- 
tinued until recently [4-7], among others Anal- 
ogous measurements under high confining pressure 
have also been reported [8—11] None of the con- 
ductivity measurements to date have studied the ef- 
fects of the oxygen fugacity or of a systematic varia- 
tion of the chemical composition These variables are 
known to have a strong influence on the conductivi- 
ties of minerals. For example, Duba et al. [1 ] and 
Duba and Ito [32] have shown that the conductivity 

TABLE 1 

Electron probe analyses {wt.%) of glasses quenched from post-run melts of tholeiite basalt (70-15, PG16), alkali olivine basalt 
(BCR2), latite (V-31) and andesite (HA) Compositions HA(2N), HA(3N), HA(3F), HA(6I') and HA(lOr) correspond to base 
composition HA with additions of 2 1% Na20, 3 2%'Na20, 3 0% FeO, 6 0% FeO and 10% FeO respectively Total Fe is arbitrarily 
given as FeO 



70-15 

PG16 

BCR2 

V-31 

HA 

HA(2N) HA(3N) HA(3F) 

HA(6H) 

HA(IOF) 

S 102 

47 1 

47 5 

53 9 

61 8 

57.9 

57 6 

57 0 

57 1 

55 4 

53 4 

T 102 

3 03 

1 54 

2 58 

0 79 

0 63 

0 63 

0 62 

0 62 

0 60- 

058 

A 1203 

15 2 

15 1 

13 3 

14 2 

19 0 

18 9 

18 6 

18 7 

18 2 

17 5 

FeO 

12 2 

14.6 

13.3 

9 1 

4.99 

4 96 

4.90 

7 83 

10.4 

13 7 

MnO 

0.21 

0.13 

0 02 

0.15 

0 09 

0.09 

0 09 

0 09 

0 09 

0 08 

MgO 

7 58 

6 91 

3 57 

0.38 

3 10 

3 08 

3 05 

3.05 

2 97 

2 86 

ao 

95 

98 

75 

3 50 

77 

76 

75 

75 

7 33 - 

7 06 

Na20 

2 59 

2.95 

3 44 

4 62 

' 4 15 

6 18 

2 18 

4 09 

3 97 

3 83’ 

K 2 O 

0.74 

0 60 

1 74 

4 46 

1 06 

1 05 

1 04 

1 04 

1 01 

0 98 

Atoms (per 100 oxygens) 










Si 

29 4 

29 6 

32 8 

36 5 

33.7 

33 3 

33,1 

33 2 

32 7 

32 1 

Ti 

1 42 

0 72 

1 IS 

0 35 

0 27 

0.27 

0 27 

0 27 

0 27 

0 26 

Al 

11 2 

11 1 

96 

99 

13 0 

12 8 

12 8 

128 

12 6 

12 4 

Fe 

6 33 

7 62 

6 75 

4 49 

2 42 

2 40 

2 38 

3 81 

5 15 

6 88 

Mn 

0 11 

0.07 

0.01 

0 08 

0 04 

0 04 

0 04 

0 04 

0.04 

0 04 

Mg 

7.04 

6 42 

3 24 

0 34 

2 69 

2 65 

2 64 

2 65 
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2 56 

Ca 

6 34 

6 57 

4 90 

2 21 

4 77 

4 71 

4 68 

4.70 

4 64 

4 55 

Na 

3 13 

3 57 

4.06 

5.29 

4.67 

6 93 

8 08 

4 61 

4 54 

4 46 

K 

0 59 

0 48 

1 35 

3 36 

0 79 

0 77 

.0 77 

0.77 

0 76 

0 75 

Si + Ti + Al 

42 0 

41 4 

43.6 

46.7 

46 9 

46 4 

46 1 

46 3 

45.6 

44 8 


19 8 

20.7 

14 9 

712 

9.9 

98 

9.7 

11 2 

12.4 

14 0 

SM* 

3.72 

4 05 

5 41 

87 

5.46 

7 70 

8.9 

5.38 

5 30 

5 21 


of olivine varies over several decades with changing 
Fe^'^'/Fe^’'’ and total Fe content We present here the 
results of an experimental investigation of the elieefs 
of oxygen fugacity, temperature and composition on 
electrical conduction in magmatic liquids 

2. Experimental technique 

The rock samples listed m Table 1 were fused, 
quenched and ground several times to provide homo- 
geneous starting materials Runs were carried out 
while the liquids were suspended in a furnace on Pt- 
Rh wire circular loops bent to approximately 3 5 mm 
radii The loop assemblies provided maximum surface 
exposure to the ambient furnace atmosphere while 
minimizing reaction with the container. Furnace 
temperature was measured with two alumina-sheathed 
Pt-Rh thermocouples periodically calibrated against 
an N B S. calibrated reference thermocouple and, the 
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melting point of gold Uncertainties in the reported 
temperatures are less than ±1®C. The furnace at- 
mosphere was controlled by mixing CO 2 and H 2 
and periodically checked against a solid electrolyte 
oxygen cell Uncertainties m the reported oxygen 
fugacities are less than ±0. 1 logic units 

For measuring resistance, the containing wire loops 
formed one electrode, and a centered vertical wire 
formed the second. Upon fusion, the sample adhered 
to the loop and assumed an approximately spherical 
shape due to surface tension. The spherical symmetry 
of the sample arrangement minimized the effect of 
volume changes when the temperature was vaned and 
facilitated the resistance to conductivity conversion 
The conversion was obtained by geometric scale 
modelmg of each sample with a liquid of known con- 
ductivity. The precision of the resistance measure- 
ments relative to each other is better than ±0.1% 
Resistance was measured with an AC phase-sensitive 
lock-in amplifier This maximized signal-to-noise ratio 
and eliminated unwanted polanzation effects Resis- 
tances were measured from 14 to 40,000 Hz and 
found to be constant within ±2% over this frequency 
range. The values reported were taken at 2000 Hz. 


3. Results 

Electrical conductivities measured at 1400°C over 
the range of O 2 fugacities 10“^ to 10"° atmosphere 
are shown in Figs. 1 and 2, The chemical composi- 
tions of the liquids are given in Table 1. Since the 
microprobe analyses yield only total Fe concentra- 
tions we have used wet chemistry to determine Fe^"*" 
m glasses of comparable compositions that were 
quenched after equilibration at various temperatures 
and O 2 fugacities. These supporting runs indicate 
that at 1400°C the range of fugacities, 10~’ to 
10“° ’ atmosphere, results in a range of Fe^'*’/ 

(Fe^"*" + Fe^"*") of approximately 0.9 to 0.2 The re- 
sults in Figs 1 and 2 clearly demonstrate that the 
conductivity is relatively insensitive to the oxidation 
state of Fe m the melt, and that m estimating the 
conductivity of liquids in the magmatic range, / q^ 
may be ignored to a first approximation The slight 
decrease in conductivity observed at relatively high 
/02 (>10“^ atmosphere) for compositions PG16 
and BCR2 is real and is related to the higher total 
Fe content of these melts 



Fig 1 Electrical conductivity of melts as a function of/Q^ 
at r = 1400°C Although m detail the oxidation state of Fe in 
the liquid depends on the composition, it is more strongly 
dependent on temperature and/Q^ As a rough guide, Fe^'*’/ 
(Fe^’*' + Fe^'*') assumes values close to 0 90, 0 75, 0 50 and 
0 20 at/oj “ 10"* and 10“ ^ respectively 

for the range of compositions represented in this figure and 
in Fig 2. Compositions are indicated m Table 1 


The temperature dependence of conductivity for 
the various melts is shown in Figs. 3, 4 and 5 In 
view of the relative insensitivity of the conductivity 
to the oxidation state of Fe, these measurements 
were performed by varying the temperature while 
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Fig 2 Electrical conductivity of andesiUc melt (HA) with 
additions of Na 2 O as a function of at T = 1400° C. Com- 

positions are indicated in Table 1 
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3 Eiectncal conductivity of melts as a function of lecip- 
I temperature. 
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Fig 5 Electrical conductivity of andesitic melt (HA) with 
additions of FeO as a function of reciprocal temperature 


T (”C) 

1500 1400 1300 1200 



54 58 62 66 

lOVT(-K)-' 


Fig 4 Electrical conductivity of andesitic melt (HA) with 
additions of Na 20 as a function of reciprocal temperature 


the C 02 /(C 02 + H 2 ) gas mixture flowing into the 
furnace was held constant at 87 6% CO 2 . This value 
of the initial composition of gas flowing into the 
furnace results in T-fQ conditions of 1500"- 
10-*", UOD^-IO-’", 1300°-10-®‘, J200°- 
10~" and IS sufficiently reducing to maintain 
Fe^V(Fe^’*' + Fe^'*') > 0 90 over the entire-range of 
temperature During the up-down temperature cycles 
(50°C increments) resistance was monitored on a 
chart recorder and stable values were obtained within 
15 minutes at each temperature Reproduability of 
the data obtained dunng temperature cyclmg is excel- 
lent, typically within 1%, and the measureiMhl and re- 
producibility errors in all cases are too small to dis- 
play in Figs. 3, 4 and 5. 

4. Discussion 

Based on the following observations, we assume 
the electrical conductivities*of these melts to be 
ionic m character. 

(a) The measured conductivities of Na-doped an- 
desitic melts of this study show marked and regular 
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increases with increasing Na concentration, a corres- 
pondence generally associated with ionic conductivity 

(b) Presnall et al .[7] have deinonstrated-that 
basalt shows a rapid decrease in conductivity during 
the liquid-solid transition, also characteristic of ionic 
conduction. 

(c) It has been shown that Faraday’s law is obeyed 
in a variety of silicate glasses (cf. [13, 14]), In fact, 
the only reported observations of electronic conduc- 
tion in glasses were for those which were fused under 


reducing conditions and which contain large amounts 
of heavy -metal ions- [14] 

It should be noted that neglect of the l/T depen- 
dence of the pre-exponential term in the conductivity 
does not alter the slopes of the log a versus IjT curves 
significantly over the temperature range of this study 
Resulting shifts in the apparent activation energies are 
less than 0 01 eV Thus linearity of log a in reciprocal 
temperature can be properly interpreted as Arrhenius 
behavior. 



1800 1700 1600 1500 1800 1700 1600 

T{«K) T(°K) 


Fig. 6 Apparent activation energies for electrical conduction as a function of temperature These values were deteimined using 
conductivities measured at successive temperatures, ne^ectmg the pre-exponential temperature dependence (sec text) 
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Inspection of Figs 3, 4 and 5 reveals that some of 
the liquids show large departures from Arrhenius tem- 
perature dependence. Variations in the apparent activa- 
tion energies for representative compositions are shown 
in Fig. 6 . Compositions 70-15, PG16, BCR2 and 
HA(IOF) show an appreciable drop in activation energy 
as temperature increases. Similar behavior was also 
found in the conductivity measurements of Presnall et 
al. [7] and Khitarov et al. [10]. These compositions 
are characterized (Table 1) by relatively high ratios 
of Melts V31, HA. HA(2N) and HA(3N), 

with much lower EM^'*'/SM'‘' ratios, all have Arrhenius 
temperature dependence Silicate liquids with a high 
ratio presumably depend more strongly 
on the divalent cations for conduction (it is assumed 
in this discussion that the large anionic polymer spe- 
cies of the silicate melt, made up mostly of S 1 O 4 and 
AIO4 sub-units, do not contribute significantly to the 
conductivity) Due to their higher ionic potential 
(charge/radius) the divalent cations polarize the oxygens 
of the surrounding polymerized units m the melt 
structure to a greater degree than the monovalent cat- 
ions. The resulting “polarization cloud” surrounding 
the cation must be transported along with the ion in 
an applied electric field This amounts to a large “ef- 
fective mass” and a low mobihty through the large 
polymerized units of the liquid structure. The diva- 
lent cations are therefore thought to be more mobile 
in the regions between the polymerized units than 
mside them. At lower temperatures the greater poly- 
merization forces the divalent cations to migrate 
around the polymer species in long, tortuous paths. 
Hence, they are subject to increased collisions and 
molecular complexing probabilities per umt distance 
traveled parallel to the applied electric field. As the 
temperature is mcreased, thermal dissociation de- 
creases the average size of the polymers and the volume 
between them increases. As a result tlie divalent ions 
can migrate along more direct paths'with a consequent 
reduction in the apparent activation energy. Thus, the 
apparent activation energy for electrical conduction 
is expected to drop with increasing temperature in 
melts that have relatively high ratios. 

The monovalent cations, on the other hand, polarize 
the neighbonng oxygens to a lesser degree and are 
more mobile through the polymerized units of the 
melt. Consequently, their path through the melt in an 
applied field will be less dependent upon the size and 


distribution of the polymerized units, and the apparent 
activation energies should show much less temperature 
dependence. This explanation is consistent with the 
data of Figs. 3 and 4 In addition, Hakin and Uhlmann 
[15] , Otto andMilberg [16] and Haven and' Verkerk 
[17] also obtained comparable activation energies 
down to much lower temperatures in alkali silicate 
glasses. 

Melts with relatively large concentrations of Fe 
and high SM^‘’'/ 2 M+ ratios (hence appreciably depen- 
dent upon divalent cations for conduction) show a 
decrease of conductivity as the Fe^‘*‘/Fe^'^ ratio de- 
creases (cf Fig. 1). Several reasons can be offered to 
explam this observation. First, a melt with more Fe^"*" 
contains fewer Fe 10 ns free to participate in conduc- 
tion. This IS simply related to the Boltzmann distnbu- 
tion of free charge carriers in the melt. Iron-oxygen 
binding energies are greater for Fe^’*' than for Fe^'*' 
Boon and Fyfe [18] and Waff et al [19] find that 
Fe^'*’ occurs only in six-fold oxygen coordination 
while Fe^'^ Mossbauer spectra show isomer shifts char- 
acteristic of tetrahedral coordination in glasses 
quenched from natural silicate liquids. The 2+ charge 
on the ferrous ion is shared among six nearest-neighbor 
oxygens. Ferric iron, on the other hand, shares a charge 
of 3-1- among either four or six oxygens and therefore 
has from 1/6 to 5/12 greater charge per Fe— 0 bond. 

A smaller fraction of feme 10 ns will be disassociated 
at any given temperature than ferrous ions Second, 
an increase of Fe^'*' ions in Fe 04 tetrahedral sites 
may reduce the effective number of univalent and di- 
valent cations that can participate in conduction 
The main polymerizing sub-units of the melt, S1O4 
can be replaced by AlO?^ or FeO^® only if M'*' or 
M^'*’ cations are available to locally offset the charge 
imbalance from substituting Al^"*" or Fe^'*' for Si'*'*' 

This charge balancing role for M'^ and M^'*' 10 ns (anal- 
ogous to their role in the feldspar structures) ties up 
a fraction of these 10 ns and effectively reduces their 
overall mobility 

Perhaps the most obvious and significant result of 
these measurements is that the total range of magmat- 
ic compositions and temperatures (approximately 
1200— I550“C) studied here accounts for less than an 
order of magnitude change in the conductivity. This 
means that the conductivity differences to be expected 
between liquids of magmatic compositional extremes 
will always be much less than those between solids 
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(e g., olivine) and liquid silicates In the melting range 
of temperatures postulated for the mantle, the partial 
melt will always be more conductive by a factor of 
2—4 decades. As discussed by Waff [2], it is probable 
that under these conditions a relatively small amount 
of partial melt can dominate the bulk conductivity 
Furthermore, a comparison of conductivities among 
the less differentiated compositions studied here (e.g., 
70-15 and PG16) and comparable compositions 
under very different conditions of confining pressure, 
and HjO and fugacities in earlier studies [7,8,10,11] 
strongly suggests that within this compositional 
grouping the temperature dependence is dominant. If 
so. It should be possible to convert geomagnetic depth- 
soundings to geothermometiy of deep-seated magma 
sources Conductivity profiles should at the very 
least be able to tell us whether we are dealing essential- 
ly with a sohd or partially melted material, i.e , 
whether or not we are in the temperature range above 
the solidus If we can make a reasonable assumption 
about the gross compositional characteristics of the 
partial melt, it may additionally be possible to estimate 
temperatures more precisely 
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ABSTRACTS 


McKAY, G. A, and D. F. Weill (1976) Petrogenesis of KREEP. Proc. Lunar Sci. 


Conf . 7th , in press. 

Abstract - Solid/liquid distribution coefficients (weight basis) have been experimen- 
tally determined for a number of trace elements for olivine, orthopyroxene, plagio- 
claj^ gnd ilmenite. Values of distribution coefficients measured at 1200'’C and f^ of 
10 ' for liquids similar in composition to the olivine-opx-plagioclase ^2 

peritectic in the pseudoternary system (Fe,Mg) SlO.-CaAl Si 0--S10- are as follows: 
olivine: Ce = 0.010±.007, Sm = 0.015±.005, Eu = 0.015±.0D5, YB = 07033±.002, 

Cr = 1.2±.l; ortho pyroxene: Ba = 0.011±.005, Cr = 5.2±.2; anorthite: Rb = 0.017±,008, 
Ba = 0.15±.03. Values measured at 1140°C and f of ^^-12. 8 for liquids similar in 
composition to high-Ti mare basalts are as ^2 follows: Ilmenite: 

Ce = 0.0061.003, Sm = 0.0101.002, Yb = 0. 0751.005. -Cr distribution coefficients 
decrease with decreasing oxygen fugaclty for both pyroxene and olivine, but are 
larger for pyroxene than for olivine at all fugacities investigated. The variation 
of with oxygen fugacity indicates that a substantial fraction of the Cr 
is divalent at lunar oxygen fugacities. 

Major and trace element partitioning and relevant phase equilibria are used 
to investigate possible parent-daughter relationships between a number of highland 
samples and highly evolved KREEP-rich materials. Out of about 80 highland samples 
tested, 33 are found to be possible parents to the KREEP-rich materials. The 
average composition of these 33 samples is very similar to that of the Low~K 
Fra Mauro basalt (LKFM) . A model is proposed in which LKFM-type material was produced 
by fractionation of large amounts of olivine and lesser amounts of plagioclase 
(and possibly other minor phases) from undifferentiated lunar material at 'U4.4 AE. 
During the same episode of igneous activity, some LKFM material differentiated to 
form a series of more evolved members of the KREEP suite, including material of 
peritectic bulk composition in the (Mg,Fe) 2 S 10 ,-CaAl 2 Sl 20 g-Si 0 „ pseudoternary 
system. In a later igneous event at 'uS.q AE, some of this peritectic material 
further differentiated, producing more highly evolved KREEP such as sample 15386. 

The model is consistent with phase equilibria, major and trace element partitioning, 
and Rb-Sr Isotopic data. 

LINDSTROM, D. J. (1976) Experimental Study of the Partitioning of the Transition 

Metals Between Clinopyroxene and Coexisting Silicate Liquids. Ph.D. thesis. 

University of Oregon. 

Abstract - This study describes the partitioning of the transition elements Sc, Ti, 

V, Cr, Mn, Fe, Co, and Ni between silicate liquids and coexisting crystalline phases 
(emphasizing clinopyroxene but also including low-Ca pyroxene, olivine, plagioclase, 
magnetite, and pseudobrookite) . 

Equilibration experiments were preformed at one atmosphere total pressure, mix- 
tures of H 2 and CO^ being used to control oxygen fugacity. Temperatures ranged from 
1110 to 1340°C. Bulk compositions were largely chosen from the system MgCaSi 0 - 
NaLlSigOg-CaAl^Sl^Og. Aliquants were doped with a few percent transition metal oxide. 
Experiments were also run on a natural basaltic composition. Quenched phases were 
analyzed with the electron microprobe. 

The crystal /liquid distribution coefficients are independent of the transition 
metal content. Analyses of pyroxenes grown from liquids with varying concentrations 
of transition metal ions presen|^indicat|^ the types of cation substitutions on the 
pyroxene lattice. Except for Mn and Ti , all species studied appear to substitute 
primarily in the Ml (Mg) site of clinopyroxene. 

Weight ratio distribution coefficients, D., exhibit considerable compositional 
dependence, narticyilarj^y between the synthetic^compositlons chosen and the basalt. 
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When equilibrium constants are calculated according to the model of Bottinga and 
Weill (1972), the compositional dependence disappears. Activities of cations in the 
liquid are approximated by their concentrations relative to the silicate melt 
network forming or network modifying components, e.g.. 




^Si 


Si Si N . + N_. + + N., + 2N_ 

Si Ti KAl NaAl CaAl^ 


N, 


Ni 




Ni 


N„. + (N^ - N„ „ ) + N,. + N„ + 

Nx ^ Ca CaAl^ Mg Fe 


These mole fractions are combined in the equilibrium constant for the reaction forming 
the pyroxene component, e.g., 


NIO^ . + CaO. , + 2SiO_ 

^ICaSl^O ^ 

^1 


(NiCaSl^O^)^^ 


liq v^iq 


^i 


Ca 
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where the activity of NiCaSi 0 is taken to be the mole fraction of Ni on the Ml site. 

Values for D. for cllnopyroxene are in the order „ 3+^„3+^„_, 2+. 2+^„ 2+^ 

2 + 3 + 2+ i+ >Sc >Co > 

Fe -vFe '^ 3 l||>Ti^ >V ‘ '' . 

When Sc enters the Ml site of ^llnopyro^ene, the charge balance is maintained 
by the simultaneous replacement of Na for Ca'^ in the M2 site or Al"^^ for Si^^ in 
a tetrahedral site. The existence of the two independent reactions precludes the use of 
a simple equilibrixim constant to describe the partitioning. 


g^romium enters pyroxene in a similar fashion, but the very limited solubility 
of Cr in the liquid presents serious experimental difficulties so that the D values 
are only approximate. 

Vanadium exists in at least two oxidation states. V"^ predominates at low 
oxygen fugaclties near the iron/wiistite buffer. In air, the observed D drops by , , 
two orders of magnitude to as low as 0.04, suggesting that V may be present. V^ 
is by n^^means ruled out, however, 

Ti enters both tetrahedral and octahedral sites in the pyroxene, octahedral 
sites b||.ng somewhat favored. 

Mn enters both the Ml and M2 sites, apparently in near equal amounts. Detailed 
modeling requires data on intracrystalline distributions not available from microprobe 
analyses. 

At the low doping levels (usually 2% Fe) used in these experiments, Fe"^"^ and Fe'^^ 
enter cj^nopyroxg^e with equal ease. 

Co and Ni behave similarly; both can be described by the equilibrium 
constants discussed above. 


A few data on low-Ca pyroxenes are reported. Observed D's for the divalent ions 
are higher than those of clinopyroxenes while D ' s for trivalent and quadrivalent 
ions are lower. 

Similar trends are apparent in the olivine data. Partitioning of the divalent ions 

is adequately described as NiO- . + 1/2 SiO_ ^ , = (NiSi„ ..0.,) , . . 

-u^ol liq 2 liq 0.5 2“^ olivine 


K 


0.5^2 

liq .liq. 1/2 


A few data on plag|.oclases show that the transition elements are rather strongly 
excluded except for Fe 


Fe-Ti oxide minerals strongly take up the transition elements. Sc and the quad- 
rivalent ions being preferentially incorporated in pseudobrookite while the other 
ions discriminate in favor of magnetite. 



WAFF, H. S. and D. F. WEILL (1975) Electrical Conductivity of Magmatic Liquids: 


Effects of Temperature, Oxygen Fugaclty and Composition. Earth Flanet ♦ Sci . Lett . , 

28, 254-260. 

Abstract - The effects of temperature, f^ and composition on the electrical con- 
ductivity of silicate liquids have been ^2 experimentally determined from 1200 
to 1550°C under a range of f cond|_^lons sufficient to change the oxidation state 
of Fe from predominantly ^2 Fe^ to Fe"^ . Oxidation of ferrous to ferric iron 
in the melt has no measurable effect on the conductivity of melts with relatively 
low ratios of divalent to univalent cations. Under strongly oxldl^fng Renditions 
a decrease of conductivity is detected in melts with high 2M /EM ratios. 

It is concluded that for purposes of estimating the conductivity of magmatic liquids, 
f^ may be ignored to a first approximation. Both univalent and divalent cation 

2 transport is Involved in electrical conduction. Melts rgfyinR heavily on divalent 
cations for conduction, i.e. melts with relatively large EM /EM ratios, show strong 
departures from Arrhenius temperature dependence with the apparent activation 
energies decreasing steadily as the temperature increases. ^Rndx^tivities dominated 
by the univalent cations, in melts with relatively small EM /EM ratios, show 
classical Arrhenius temperature dependence. These observations are discussed in terms 
of the general characteristics of the melt structure. 

Compositional variations within the magmatic range account for much less than 
an order of magnitude variation in electrical conductivity at a fixed temperature. 
This observation, combined with previous measurements of the conductivity of olivine 
(A.Duba, H.C. Heard and R. Schock, 1974)make it possible to state with reasonable 
confidence that melts occurring within the mantle will be more conductive by 3-4 
orders of magnitude than their refractory residues. Potential applications to 
geothermometry are discussed. 


LEEUAN* W, P. (X974) Experimental Determination of Partitioning of Divalent 

Cations Between Olivine and Basaltic Liquid. Ph.D. thesis. University 

of Oregon. 

- The distribution of divalent Nl, Mg, Co, Fe, Mn, and Ca between 
olivine and natural basaltic liquid was investigated between 1070®C and 
1400 C. Relative magnitudes of the distribution coefficients decrease in 
the order listed above. Distribution coefficients defined as the ratio of 
weight fractions (D) or the theoretical equilibrium constant (K) are strongly 
temperature dependent for all of these elements. Compound distribution 
^efficients (e.g., are less temperature dependent. Linear regress- 

Iot equations. In D (or K) « A/X + B, are given for each element. 

These equations nay be used as geo thermometers to obtain olivine crystalliz- 
ation temperatures. Applications to experimental runs and natural volcanic 
rocks suggest that under optimum conditions temperatures can be estimated 
with an accuracy of better than ±30", Simple D values are apparently not 
very sensitive to differences in bulk composition within the range of 
terrestrial basalt. Temperature calibrations of D may be useful in construct- 
ing detailed crystallization or melting models involving olivine. 



WEILL, D, F. and G. A. McKAY (1975) The partitioning of Mg, Fe, Sr, Ce, Sm,' Eu, 


and Yb in Lunar Igneous Systems and a Possible Origin of KREEP by 
Equilibrium Partial Melting. Proc . Lunar Sci. Conf . 6th , 1143-1158. 


Abstract - The solid/liquid distribution of Mg, Fe, Sr, Ce, Sm, Eu, and Yb have 
been determined over crystallization paths in the system (Mg,Fe) 2 Si 0 ^-GaAl 2 Si 20 g- 
SiO^. Varying degrees (1-5 wt.%) of partial melting of an olivine + ^orthopyroxene 
+ plagioclase rock yield liquids with large-ion llthophlle (LIL) and major-element 
concentrations strikingly similar to those of some Apollo 15 and 16 KREEP rocks. 

The parent assemblages from which these KREEP melts can be derived have a bulk 
composition similar to that previously proposed for the whole moon. Partial 
melts derived from a source with the same mineralogy but with a chemical 
composition corresponding to the average lunar highlands do not simultaneously 
satisfy the major-element and LIL concentration characteristics of KREEP. 



«STU. , D. F., MCKAY, G. A.» KR1DELSAU6H, S. J. and M. 0»JIZECK (1974) 
Modaling the Evolution of Sm and Eu Abundaneea During Lunar Igneous 
Differentiation. Geochlm . Cogmochlm , Acta.S uppl. 5, vol. 2, 1337- 
1352. 

Abstract ~ Models are presented for the evolution of Eu and Sa abundances 
during lunar igneous processes. The effect of probable variations in 
lunar teaperature and oxygen fugacity, mlDeral->llquld distribution coeffi' 
dents » and the crystallization or melting progression are considered la 
the model calculations. Changes In the proportions of crystallizing phases 
strongly influence the evolution of trace element abundances during frac- 
tional crystallizatiODy and models must include realistic estimates of the 
maj^or phase equilibria during crystallization. 

The results are applied to evaluate the possibility of generating 
KREEP-rlch materials by lunar igneous processes. Fractional crystalliza- 
tion of a primordial liquid siiell is a very unlikely origin because the 
characteristic Sm and Sm/£u abundances of KREEF are achieved only after a 
very high degree of fractionation involving separation of mostly ferro- 
magnesian minerals. The low Fe/Mg found in KREEF is incompatible with 
such a history. All fractional crystalllzatlbn sequences involving 
appreciable separation of plagloclase produce negative Eu anomalies that 
are much too large. 

Partial melting of ollvlne-orthopyroxene-cllnopyzoxene-plagloclase 
assemblages can produce KREEF abundances of Sm and Sm/Eu. The most 
probable parent assemblages are (1). dominantly olivine plus orthopyroxene 
(<20Z plagloclase and a few percent cllnopyroxene) with bulk Sm and Sm/Eu 
concentrations comparable to those estimated for the whole moon; i.e. 
enriched approxtmatsly tenfold or less relative to chondrites, -with no Eu 
anomaly, or (2) a rock with more plagloclase, eiuriched in both ‘Sm and Eu 
(no negative Eu anomaly) relative to the whole moon. An origin for either 
of these parental bulk compositions via fractional crystallization implies 
a low concentration of plagloclase constituent relative to ferrcxsagneslan 
constituents In the moon's primordial liquid shell. 

Although the range of Sm and Sm/Eu ^undances may be related to vary- 
ing degrees of partial melting of these parental assemblages, it is also 
possible that the variation may simply reflect different amounts of olivine 
plus ortbopyroxene relative to plagloclase plus cllnopyroxene in the parent 
material. 

DRAKE, H. J. (1975) The oxidation state of europium as an indicator of 

oaygen fugacity. Geochlm . Cosmochlm . Acta. 39, 55-64. 

Abstract - The distribution of Eu between plagloclase feldspar and magmatic 
liquid has been determined experimentally for basaltic and andesitic systems 
as a function of temperature and oxygen fugacity at one affioap^re total 
pressure. Using the approach of Phllpotts, the ratios Eu'"/£u^ in plagio- 
clase and coexisting magmatic liquid have been calculated. These ratios 



appear co be almply related to oxygen fugaclty for the bulk compositions 
studied here. Using published trace element distribution data for natural 
rocks » oxygen fugacltles may be calculated from these experimental results. 

For terrestrial basalts, calculated oxygen fugacltles average 10“^ with 
little dispersion from this value* Andesites average lO"^*^ with consider- 
able dispersion, while dacites and rhyodacites average also with 

considerable dispersion. Oxygen fugacltles for lunar ferrobasalts cluster 
tightly around 10“' Data on achondrltic meteorites are limited, but 
calculations Indicate oxygen fugacltles of two-to~£ive orders of magnitude 
lower than lunar ferrobasalts. 

DUHCAN, A. R., GSIRVE, R. A. and D. F. HEILL (1975) The life and times of 
Big Bertha: lunar polymlct breccia 14321. Geochim . Cosmochim . Acta , 

39. 265-273. 

Abstract - The lithlc units of polymlct breccia 14321* Bertha) have been 
grouped according to composition, texture, degree of metamorphism, and addi- 
tional criteria based on a systematic study of the interrelationships of all 
clast-matrix pairs. From this information it has been possible to reconstruct 
the assembly emd metamorphlc history of this breccia. The earliest formed 
fragmental component of 14321 (mlcrobreccia-1) Is dominated by KREEP-rich 
fiorlte, extruded and subsequently brecclated and llthlfied in an ejecta blan- 
ket at approximately 1000^ C In the general region of Mare Imbrium after the 
Serenltatls impact but prior to the Imbrium Impact. This early mlcrobreccla 
component and lesser amounts of mare-type basalt, mlcrogranite, rhyolite 
glass, anorthosite and olivine mlcrobreccla were assembled at the Apollo 14 
site as part of the Fra Mauro ejecta blanket from the Imbrium impact. The 
resulting mlcrobreccia-3 incorporates all the lithlc types above and accre- 
tionary lapilli structures (mlcrobreccia-2) In a dark matrix annealed at 
approximately 700°C. A later Impact on the Fra Mauro excavated and mutually 
abraded mlcrobreccia-3 and a local, 14321-type, basalt which were assembled 
Into polymlct breccia 14321. Final placement of 14321 at its sampling 
location was accomplished during the minor Cone Crater impact event* 

GRIEVE, R. , HCKAY, G. , SMITH, H. and D. F. WEILL' (1975) Lunar polymlct 
breccia 14321: a petrographic study. Geochim . Cosmochim. Acta , 39, 

229-245. 

Abstract - The lithlc units of polymlct breccia 14321 (Big Bertha) have been 
grouped according to composition, texture, degree of metamorphism, and 'addi- 
tional criteria based on a systematic study of the interrelationships of all 
clastHnatrix pairs. From this Information It has been possible to reconstruct 
the assembly and metamorphlc history of this breccia. The earliest formed 
fragmental component of 14321 (mlcrobreccla-1) Is dominated by KREEF-rich 
norite, extruded and subsequently brecclated and llthlfied in an ejecta blan- 
ket at approximately 1000° C In the general region of Mare Imbrium after the 



Serenltacls impact but prior to the Imbriiim impact. This early microbreccla 
component and lesser amounts of mare-type basalt, microgranite, rhyolite 
glass, anorthosite and olivine microbreccia were assembled at the Apollo 14 
site as part o£ the Fra Mauro ejecta blanket from the Imbrlum impact* The 
resulting mlcrobreccia-3 incorporates all the llthic types above and accre- 
tionary lapilll structures (microbreccla-2) in a dark matrix annealed at 
approximately 700°C. A later Impact on the Fra Mauro excavated and mutually 
abraded microbreccia-3 and a local, 14321-type, basalt which were assembled 
into polymlct breccia 14321. Final placement of 14321 at its sampling loca- 
tion was accomplished during the minor Cone Crater impact event. 

DRAKE, M. J. and D* F. WEILL (1975) Partition of Sr, Ba, Ca, Y, and 

Eu and other REE between plagloclaae feldspar and magmatic liquid: an 

experimental study. Geochlm . Cosmochim . Acta , in press* 

Abstract - Plagioclase feldspar /magmatic liquid partition coefficients for Sr, 
Ba, Ca, Y, Eu^*^, £u^+ and other REE have been'deteimlned experimentally at 
one atmosphere total pressure in the temperature range 1150-1400^0. Natural 
and synthetic melts representative of basaltic and andesitic bulk compositions 
were used, crystallizing plagioclase feldspar in the composition raage An $5 
to Ang^. Partition coefficients for Sr are greater than unity at all geologi- 
cally reasonable temperatures, and for ‘Ba are less than unity above approxi- 
mately 1060°C. Both* are strongly dependent upon temperature. Partition 
coefficients for tHe trlvalent REE are relatively Insensitive to temperature. 
At fixed temperature they decrease monotonlcally from La to Lu. The parti- 
tion of £u is a strong function of oxygen fugaclty. Under extreme reducing 
conditions 0^^ approaches the value of Ug^* 

ORUTZECK, M. , KRIDELBAUGU, S. and D.' F. WEILL (1974) The distribution of 
Sr and REE between dlopslde and silicate liquid. Geophys . Res . 

Letters , 1, 273-275. 

Abstract - Distribution coefficients (weight concentration in the solid 
divided ly concentration in the liquid) have been detensined experimentally 
In diops ide-liquid pairs. Experiments were carried out In air at 1265*C in 
the system NaAlSi308— CaAl 2 Si 20 g— CaMgSi 20 g at the two bulk compositions: 

5-30-65 and 32-11-57 (weight percentage of components in the order given 
above) . Small differences are observed in the distribution coefficients 
between the two compositions. Mean values of the coefficients In the first 
system are; Sr=.078, La-.069, Ce-.098, Nd«,21, Sm-.26, Eu-.31, Gd«.30, 

Er».30 and Lu».28. At low oxygen fugaclties of 10”®»5 and 
atmospheres the distribution of Eu between dlopslde and liquid develops a 
negative Eu-anomaly‘‘ with coefficients of .23 and .12 respectively. 



SMITH » H. D. <1973) An Experimental Study of the Diffusion of Na» K, and 

Rb in Hagaatic Silicate Liquids. Ph.D. thesis. University of Oregon. 

Abstract - In this study an experimental method was developed for the 
determination of diffusion coefficients in multicomponent silicate melts. 
This method was then applied to the study of the diffusion of Na*^, K+, and 
Rb*^ in three silicate melts which were modeled composltionally after a 
basalt, an andesite, and a rhyolite. 

Diffusion in a multicomponent system is described generally by a 
"diffusion matrix" [see Onsager (1945), Cooper (1965)]. However, subject 
to certain limitations, a single coefficient can be used to describe the 
diffusion of a single component through a liquid of given composition. 

This coefficient Is the Effective Binary Diffusion Coefficient, ££DC, 
described by Cooper (1968). 

In this study the EBDCs for Na***, id", and Kb+ were determined by setting 
op diffusion couples in which concentration gradients were present only 
for the alkalies. Cylinders of melt with the same base melt composition 
but containing different alkalies were joined end-to-end and diffusion 
allowed to take place. After quenching, each couple was sectioned length- 
wise and an electron microprobe was employed to determine the concentration 
profiles that developed for each alkali. From these concentration profiles 
it was possible to determine the EBDC at regular intervals along the profile 
by graphical application of the Boltzmann-Matano analysis. The D.s (EBDC) 
so determined cover a wide range of temperature (1128* C to 1535*C7 and 
compositional variables. 

The findings of this experimental study can be summarized as follows: 

1) Stroeg electrical coupling between ionic fluxes is indicated by the near 
equality of opposed alkali fluxes in most couples. A measurable amount of 
"uphill diffusion" of calcium in some couples of basaltic and andesitic 
(calcium-rich) compositions also demonstrates the ability of other relatively 
mobile ions to participate in ionic exchange dlffuslbn within silicate melt 
structures. 2) D^ varies with position along the concentration profile of 
1 (conveniently represented by the ratio w, 5 C.-Cj/Cj-Ci in which C, is the 
concentration of 1 and and C 2 are the concentrations at the ends of the 
effectively infinite couple). For example, there is a factor of 20 variation 
for Djy^ and Dj. between w « 0 and u ■ 1 in Rb-Na-rhyollte couples. In con- 
..trast, and^„ exhibit very little variation across the Rb-K-basalt 
couples. These two examples represent f:he extreme cases found in this study. 
Comparisons of the diffusion coefficients of the alkalies for the purpose of 
evaluating the effects of temperatures or matrix composition must be made 
at comparable u values (u 'f* in this study) . 3) The D , (EBDC) values are 

Arrhenius functions within the range of temperature used in this study. The 
apparent activation energies, E., vary from 26’ Kcal/mole for D^ in K-Na- 
xhyollte to 48 Kcal/mole for D . in Rb-K-basalt systems. 4) E^^varies 
systematically with the liquid matrix composition. A linear rilationshlp 
was found between and the percentage of silicate network forming (HHF) 
units in the liquid for each ion pair (Rb-Na, K-Na, Rb-K). £, decreases with 
increasing NWF percentage. 5) Analysis of £ and B* via absolute rate theory 
suggests that rhyolite compositions form open structured melts while basalt 
compositions form relatively compact-structured melts. The concurrent 
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variation of E. and (the “Compensation" effect) Is In full agreement vrlth 
that shown by alf fusion data In other silicate systems [see Winchell (1969)]. 
6) fiaealt^andesite nixed couples demonstrate that the D.s for the alkalies 
are much larger (about 30 times) when exchange between different alkalies 
can occur. When only one alkali is present in a couple, the D^s for all the 
elements are about the same magnitude. 

In the final section, three cases of diffusion in natural systems are 
analysed, two of lunar and one of terrestrial origin, 

WEILL, D. F. and H, J. DRAKE (1973) Europium Anomaly in Flagioclase: Experi- 
mental Results and Semiquantitative Model, Science, 180, 1059-1060. 

Abstract - The partition of europium between plagioclase feldspar and magmatic 
liquid is considered in terms of the distribution coefficients for divalent 
and trivalent europium. A model equation is derived giving the europium 
anomaly in plagioclase as a function of temperature and oxygen fugaclty. The 
model explains europium anomalies in plagioclase synthesized under controlled 
laboratory conditions as well as the variations of the anomaly observed in 
natural terrestrial and extraterrestrial igneous rocks. 



